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Preface 


This  final  technical  report  describes  a  three-year  research  program  to  investigate 
computer  generated  holograms  written  by  electron  beam  lithography.  This  work  was 
performed  at  Honeywell  Corporate  Technology  Center  under  sponsorship  of  the  Air 
Force  Office  of  Scientific  Research,  Contract  No.  F496208O-C-OO29.  The  period  of  con¬ 
tract  performance  was  from  1  January  1980  to  28  February  1983.  The  contract 
monitor  has  been  Dr.  John  A.  Neff,  (former)  program  manager  of  Electronics  and 
Material  Science,  AFOSR/NE.  The  principal  investigator  is  Steven  M.  Arnold. 


Abstract 


This  final  report  describes  a  three-year  research  program  to  investigate  computer 
generated  holograms  produced  by  electron  beam  lithography.  Because  they  are  light 
weight  and  can  create  wavefronts  of  arbitrary  complexity,  computer  generated  holo¬ 
grams  are  attractive  for  applications  such  as  optical  data  processing  and  optical  tes¬ 
ting.  To  be  practical,  these  elements  must  have  large  space-bandwidth  products,  and 
have  the  qualities  of  high  diffraction  efficiency,  low  scattering,  low  aberration  and 
low  cost,  as  compared  to  conventional  optics.  We  have  developed  generalized  encod¬ 
ing  algorithms  and  fabrication  techniques  for  producing  e-beam  computer  generated 
holograms  having  submicron  feature  sizes,  distortion-free  resolution  of  better  than 
0.4  micron,  and  space-bandwidth  products  in  excess  of  107.  E-beam  lithography  is  su¬ 
perior  to  optical  plotting  for  writing  computer  generated  holograms  and,  with 
further  refinement,  could  be  extended  to  produce  holograms  having  space-bandwidth 
products  as  large  as  1011,  comparing  favorably  with  interferometrically  recorded 
holograms. 

We  have  demonstrated  the  application  of  e-beam  computer  generated  holograms  to 
the  testing  of  asymmetric,  aspheric  optics  having  up  to  123  waves  of  aspheric  aber¬ 
ration.  We  have  also  developed  an  optical  vector-matrix  multiplication  scheme  that 
utilizes  e-beam  computer  generated  holograms  to  achieve  high  numerical  accuracy 
and  large  numerical  range;  an  accuracy  of  20  dB  and  a  range  of  37  dB  have  been 
demonstrated. 

Techniques  for  increasing  diffraction  efficiency  beyond  the  10.1-percent  limit  of 
binary  absorption  holograms  have  been  investigated.  Ion  milling  has  been  used  to 
produce  40-percent  efficient,  environmentally  durable,  all-glass  holograms.  Tech¬ 
niques  for  producing  blazed,  reflection  holograms  by  anisotropic  etching  of  silicon 
have  also  been  explored.  The  further  development  of  high-efficiency,  large  space- 
bandwidth  product,  e-beam  computer  generated  holograms  is  expected  to  have  a 
marked  impact  on  the  future  of  holographic  optics. 
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Section  1 
Introduction 


Holographic  optical  elements  can  play  an  important  role  in  many  areas  of  optical  en¬ 
gineering.  Because  they  work  by  diffraction,  they  offer  many  advantages  over  con¬ 
ventional,  refractive  optical  elements,  most  notably:  compactness,  light  weight,  and 
low  replication  cost.  Against  these  advantages  must  be  weighed  their  high  dis¬ 
persion  (chromatic  aberration)  and  limited  diffraction  efficiency. 

Many  research  groups  have  demonstrated  holographic  optical  elements  in  practical 
optical  systems.  Holographic  lenses  have  been  employed  in  interferometers,  wherein 
they  provide  a  reference  wavefront  for  comparing  optical  surfaces  [1-1].  They  have 
also  been  used  in  laser  machining  to  shape  the  laser  [1-2].  In  1979,  Fienup  and 
Leonard  [1-3]  demonstrated  that  optical  processors  utilizing  holographic  optics  were 
capable  of  processing  images  of  large  space-bandwidth  product.  More  recently, 
Horner  and  Ludman  [1-4]  have  shown  the  use  of  holographic  optical  elements  for  de¬ 
multiplexing  in  fiber  optic  systems.  To  be  more  compatible  with  existing  optical  sys¬ 
tems,  these  elements  must  have  the  qualities  of  low  aberration,  high  diffraction  effi¬ 
ciency  over  the  desired  field  of  view,  low  scattering,  and  low  cost. 

With  ordinary  interferometric  recording  of  holographic  optics,  we  are  limited  to  the 
complexity  and  accuracy  of  available  recording  wavefronts,  and  also  by  the  avail¬ 
ability  of  suitable  recording  materials.  If  the  recording  and  design  wavelengths  dif¬ 
fer,  creation  of  these  recording  wavefronts  can  become  exceedingly  difficult. 

With  computer  generated  holography,  recording  wavefronts  are  no  longer  necessary. 
The  desired  interference  pattern  is  instead  calculated  and  then  plotted  using  a  suit¬ 
able  graphic  device.  The  advance  of  computer  generated  holography  has  opened  a 
new  frontier  in  holographic  optics.  It  is  by  far  the  best  technique  for  making  optical 
elements  to  form  wavefronts  of  arbitrary  complexity. 

Techniques  for  making  binary  computer  generated  holograms  are  well  developed. 
Their  potential  in  optical  engineering  applications  is  ever  increasing.  Computer  gen¬ 
erated  holograms  have  been  used  in  such  areas  as  optical  data  processing,  optical 
testing,  optical  memories,  laser  beam  scanning,  and  3-D  image  display.  An  article  by 
W.H.  Lee  [1-5]  is  an  excellent  source  of  information  on  both  the  techniques  of  mak¬ 
ing  computer  generated  holograms  and  their  applications.  With  the  development  of 
improved  fabrication  techniques,  more  applications  of  computer  generated  holo¬ 
grams  are  anticipated  in  thufature. 

The  computer  generation  of  a  hologram  involves  the  following  three  major  steps: 

•  Calculation  of  the  complex  wavefront  at  the  hologram  plane  originating  from  a 
given  wavefront  at  the  object  plane.  For  holographic  lenses,  this  wavefront  is 
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The  typical  procedure  for  fabricating  a  computer  generated  hologram  is  to  have  the 
digitized  interference  pattern,  which  has  been  calculated  and  encoded  by  computer, 
drawn  to  a  large  scale  by  a  computer-driven  plotter.  The  drawing  is  then  reduced 
photographically  onto  high-resolution  film  to  the  desired  final  size.  There  are  sever¬ 
al  disadvantages  and  inherent  limitations  with  this  procedure:  1)  errors  are  in¬ 
troduced  in  the  plotting  and  photo-reduction  processes;  2)  optical  plotting  devices  are 
limited  in  spatial  resolution  and  space-bandwidth  product,  typically  to  106  pixels; 
and  3)  the  turnaround  time  for  the  indirect  procedure,  from  plotting  to  photo-reduc¬ 
tion  to  film  development,  can  be  days,  which  is  impractical  for  industrial  applica¬ 
tions. 

Electron  beam  lithography  overcomes  all  of  the  above  disadvantages  and  limita¬ 
tions.  E-beam  writing  of  computer  generated  holograms  is  depicted  in  Figures  1-2 
and  1-3.  Distortion  errors  are  significantly  reduced  because  of  direct  writing  and  ex¬ 
cellent  spatial  resolution.  Figure  1-4  shows  the  submicron  feature  capability  of  holo¬ 
grams  written  by  e-beam  lithography.  Holograms  with  linewidths  as  small  as 


Figure  1-2.  E-Beam  Direct-Writing  of  a  Computer  Generated  Hologram.  The 

fringe  pattern  is  created  by  a  combination  of  beam  deflections  and 
workstage  translations.  Both  are  under  control  of  the  e-beam 
computer,  which  may  also  function  to  generate  the  fringe  pattern 
data. 
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^E-BEAM  RESIST  (3000-10.000A) 
CHROME  ( -  500A) 

GLASS  OR  QUARTZ  SUBSTRATE 


EXPOSE  HOLOGRAM  PATTERN 
WITH  COMPUTER  CONTROLLED 
E-BEAM 


DEVELOP  RESIST 


ETCH  THROUGH  CHROME 


Fig  ire  1-3.  Fabrication  Process  for  Binary,  Chrome-on-Glass  Type,  Computer 
Generated  Holograms.  These  holograms  can  have  diffraction 
efficiencies  of  10.1  percent. 


IOOoA  have  been  demonstrated  [1-6].  By  comparison,  the  most  advanced  optical 
hologram  writer  can  produce  linewidths  no  narrower  than  5  microns  [1-7].  For  e- 
beam  computer  generated  holograms,  the  achievable  number  of  pixels  ( >  1010)  can 
approach  that  of  interferometrically  recorded  off-axis  holograms  because  of  the  sub¬ 
micron  resolution  and  the  capability  of  having  many  small  scan  fields  (typically  1 
mm  x  1  mm)  stitched  together  by  interferometrically  controlled  translation  of  the 
workpiece.  An  advanced  e-beam  lithography  system  can  easily  produce  a  hologram 
directly  onto  a  3-inch  substrate  in  45  minutes,  offering  a  much  shorter  turnaround 
time  for  practical  applications  when  compared  to  indirect  plotting. 
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Figure  1-4.  Grating  Pattern  in  Honey  Wed-Proprietary  Electron  Resist  EP25,  as 
Produced  by  Honeywell's  Cambridge  EBMF  E-Beam  Lithography 
System.  The  grating  period  is  1.5  microns. 


RESEARCH  OBJECTIVES 

The  herein  reported  work  under  AFOSR  Contract  No.  F49620-80-C-0029  initiated  e- 
beam  computer  generated  hologram  research  at  Honeywell  beginning  on  1  January 
1980. 

Our  overall  long-range  objective  in  this  research  has  been  to  investigate  and  develop 
unique  e-beam  computer  generated  holograms,  and  to  advance  the  fundamental 
knowledge  of  computer  generated  hologram  technology  for  applications  in  optical  re- 
'  search  of  current  and  future  interest. 

Our  short-term  objectives  for  the  1980  research  program  were: 

•  To  evaluate  the  quality  of  e-beam  writing 

•  To  demonstrate  the  feasibility  of  using  e-beam  lithography  to  write  high-per¬ 
formance  computer  generated  holograms  for  use  in  aspheric  testing. 


Our  initial  efforts  concentrated  on  evaluating  the  quality  of  e-beam  writing  using  an 
existing  Honeywell  e-beam  system  and  a  plane  wave  interference  technique.  These 
studies  not  only  advanced  our  understanding  of  making  computer  generated  holo¬ 
grams  by  e-beam  direct-writing,  but  also  enabled  us  to  improve  our  e-beam  writing 
accuracy  by  software  compensation.  At  the  same  time,  we  performed  a  distortion 
analysis  of  a  state-of-the-art  optical  plotting  device.  The  results  of  this  analysis  con¬ 
firmed  the  superiority  of  e-beam  lithography  for  writing  computer  generated  holo¬ 
grams,  and  quantified  the  limitiations  of  indirect  optical  plotting  techniques. 

During  the  program,  optimized  software  was  developed  to  efficiently  encode  aspheric 
wavefronts  of  arbitrary  complexity.  We  produced  a  computer  generated  hologram  of 
an  symmetric  aspheric  wavefront,  designed  to  be  tested  using  the  Honeywell/Tropel 
holographic  and  shearing  interferometer,  thus  demonstrating  the  feasibility  of  writ¬ 
ing  complex  holograms  by  e-beam  lithography.  This  technology  has  since  been  used 
to  test  aspheric  diamond-turned  optics  at  Honeywell’s  Electro-Optics  Division. 

As  a  continuation  of  our  efforts  to  develop  unique  holographic  optical  elements  using 
e-beam  lithography,  we  proposed  for  1981  to  investigate  and  develop  a  new  kind  of 
hologram,  the  partitioned  computer  generated  hologram  (PCGH),  which  is  useful  in 
optical  computing. 

Our  short-term  objectives  for  the  1981  research  program  were: 

•  To  analyze  and  develop  capabilities  for  fabricating  PCGHs  using  the  unique 
capabilities  of  e-beam  direct-writing. 

•  To  demonstrate  PCGHs  applicable  to  optical  computing,  with  maximum  numer¬ 
ical  range  >  40  dB  and  maximum  accuracy  >  40  dB. 

The  PCGH  differs  from  conventional  computer  generated  holograms  because  it  is 
spatially  partitioned,  fully  uses  available  light  by  means  of  a  space-variant  illumina¬ 
tion  profile  (such  as  a  Gaussian  distribution),  and  allows  trimming  for  amplitude 
control.  The  unique  capability  of  e-beam  writing  is  especially  useful  to  the  PCGH. 

During  1981,  we  developed  an  approach  to  the  design  of  PCGHs,  calculated 
achievable  values  of  numerical  range  and  accuracy,  and  implemented  partitioning 
algorithms  to  produce  PCGHs  by  e-beam  lithography.  We  designed  and  fabricated 
PCGHs  to  perform  the  optical  equivalent  of  a  1  x  10  matrix  operation.  In  separate 
experiments,  we  encoded  matrices  with  20  dB  of  numerical  accuracy  and  37  dB  of 
numerical  range.  These  results  were  limited  primarily  by  the  stability  of  the  laser 
light  source. 

Computer  generated  holograms  are  by  far  the  best  method  for  making  optical  ele¬ 
ments  to  form  wavefronts  of  arbitrary  complexity,  but  their  use  has  not  been  wide¬ 
spread.  The  reasons  for  this  are  several.  They  include  limited  space-bandwidth  prod¬ 
uct,  large  computational  effort,  and  low  diffraction  efficiency.  We  have  shown  that  e- 
beam  lithography  can  largely  overcome  the  space-bandwidth  limitation,  and  we 
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have  demonstrated  versatile  and  efficient  encoding  algorithms.  As  a  continuation  of 
our  efforts  to  develop  useful  computer  generated  holograms  using  e-beam  litho¬ 
graphy,  we  proposed  for  1982  to  investigate  and  develop  techniques  for  fabricating 
high-efficiency  computer  generated  holograms. 

Our  short-term  objectives  for  the  1982  proposed  research  program  were: 

•  To  analyze  and  develop  techniques  for  fabricating  e-beam  computer  generated 
holograms  with  diffraction  efficiencies  approaching  100  percent. 

•  To  demonstrate  high-efficiency,  large  numerical  aperture,  computer  generated 
holograms  capable  of  bringing  a  collimated  He-Ne  laser  beam  to  a  focus. 

-The  diffraction  efficiency  of  our  chrome-on-glass  computer  generated  holograms  is 
limited  to  a  theoretical  maximum  of  10.1  percent.  We  have  investigated  two  ap¬ 
proaches  for  achieving  higher  efficiency.  The  first  approach  involves  ion  milling  of 
the  hologram  pattern  into  the  glass  to  produce  an  all-glass,  surface-relief  type  grat¬ 
ing.  Such  a  hologram  is  very  environmentally  durable  and  has  a  high  damage 
threshold.  We  designed  and  fabricated  several  ion-milled  holograms  with  diffraction 
efficiencies  of  40  percent.  While  we  had  hoped  to  achieve  80-percent  efficiency  as 
grating  periods  approached  a  micron,  difficulties  in  controlling  linewidth  in  the 
thick  resists  prevented  such  gains. 

Our  second  approach  to  high-efficiency  computer  generated  holograms  utilizes 
anisotropic  etching  to  produce  blazed  grating  patterns  in  silicon.  We  designed  and 
fabricated  a  number  of  such  blazed,  crystalline  holograms  to  either  focus  a  plane 
wave  or  collimate  a  spherical  wave.  These  holograms  rely  on  sharply  peaked  groove 
profiles  for  their  diffraction  efficiency.  Problems  with  undercutting  of  the  etch  mask 
adversely  affected  the  groove  profiles,  preventing  diffraction  efficiencies  in  excess  of 
about  20  percent. 

Despite  our  limited  success  in  fabricating  them,  we  feel  that  e-beam  produced  sur¬ 
face-relief  holograms  have  considerable  potential. 

The  following  sections  report  the  program  work  in  detail.  Section  2  discusses  the  use 
of  e-beam  lithography  in  the  writing  of  computer  generated  holograms,  and  presents 
our  hologram  encoding  methods.  Section  3  reports  our  application  of  e-beam  com¬ 
puter  generated  holograms  to  problems  of  aspheric  optical  testing  and  optical  vector- 
matrix  multiplication.  Section  4  reports  our  efforts  to  develop  fabrication  techniques 
for  increased  diffraction  efficiency.  Section  5  summarizes  the  program  accomplish¬ 
ments. 
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Section  2 

Electron  Beam  Fabrication  of  Computer  Generated  Holograms 


Two  essential  requirements  must  be  met  to  achieve  a  useful  computer  generated  hol¬ 
ogram.  First,  an  adequate  graphic  device  is  needed  for  translating  the  conceptual 
and  digital  representation  of  the  hologram  into  a  recording  medium.  The  graphic  de¬ 
vice  must  be  capable  of  plotting  hundreds  of  lines  constructed  from  thousands  of  line 
segments  to  an  accuracy  of  a  fraction  of  the  average  line  spacing.  Commercially 
available  optical  recorders  can  provide  only  about  2000  x  2000  distortion-free  pixels 
across  the  hologram.  This  is  far  below  the  1010  or  more  pixels  of  a  typical  conven¬ 
tional  hologram. 

Second,  adequate  software  is  needed  to  drive  the  graphic  device.  The  software  will 
have  to  execute  plotting  procedures  millions  of  times  and  must  be  efficient  and  flex¬ 
ible  to  economize  the  total  plotting  time.  Preprocessing  of  data  can  reduce  plotting 
time,  but  a  compact  data  format  is  needed  if  data  storage  requirements  are  to  be 
kept  manageable. 

Our  overall  objective  in  this  program  has  been  to  extend  the  use  of  e-beam  litho¬ 
graphy  to  the  making  of  computer  generated  holograms.  In  e-beam  lithography,  a 
highly  collimated  beam  of  electrons,  under  computer  control,  is  used  to  expose  a  pat¬ 
tern  in  electron  resist.  The  area  over  which  the  beam  can  be  deflected  with  precision 
is  limited  (1  mm  x  1  mm  is  typical)  so  that  large  patterns  are  made  by  translating 
the  workpiece  and  stitching  small  pattern  segments  together. 

Two  approaches  to  the  building  of  large  patterns  have  been  developed.  One  is  where 
the  workpiece  is  moved,  brought  to  rest  and  the  pattern  segment  exposed,  usually  in 
a  vector-scan  technique.  The  two  lithography  systems  used  in  the  present  work  were 
of  this  type.  Early  phases  of  the  program  used  a  Honeywell-developed  electron  beam 
exposure  system  (Honeywell  EBES).  Tests  showed  this  instrument  to  be  far  superior 
to  a  state-of-the-art  optical  plotting  device.  In  late  1980,  Honeywell  obtained  an 
EBMF  e-beam  lithography  system  from  Cambridge  Instruments.  Because  of  its  supe¬ 
rior  resolution,  this  system  soon  supplanted  the  Honeywell  EBES  for  hologram 
drawing. 

A  second  approach  to  pattern  composition  is  to  keep  the  workpiece  in  continual  mo¬ 
tion  while  exposing  a  small  stripe  of  the  pattern  in  a  raster-scan  technique.  Honey¬ 
well  has  recently  acquired  a  system  of  this  type,  based  on  Bell  Lab’s  EBES  design. 

This  section  discusses  significant  aspects  of  the  Honeywell  EBES  and  Cambridge 
EBMF  e-beam  lithography  systems  and  compares  e-beam  lithography  to  optical  plot¬ 
ting.  Our  hologram  encoding  methods  are  presented.  The  section  concludes  with  a 
discussion  of  the  capabilities  and  limitations  of  the  Cambridge  EBMF  as  a  hologram 
writing  tool. 

I 
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HONEYWELL  ELECTRON  BEAM  EXPOSURE  SYSTEM  (HONEYWELL  EBES) 

* 

Honeywell  initiated  an  e-beam  lithography  program  in  1974  that  included  both  sys¬ 
tem  development  and  electron  resist  research.  Figure  2-1  is  a  photograph  of  the 
Honeywell  EBES  system.  In  the  following  paragraphs,  pattern  distortion  of  the 
Honeywell  EBES  is  discussed,  comparison  is  made  to  a  state-of-the-art  optical  re¬ 
corder,  and  our  hologram  encoding  method  is  presented. 

Distortion  Analysis 

To  demonstrate  the  feasibility  of  e-beam  computer  generated  holograms,  a  prere¬ 
quisite  is  to  determine  the  distortion-free  resolving  power  of  the  e-beam  system. 
Then  the  system  can  be  used  to  make  holograms  of  predictable  optical  figure.  Pat¬ 
tern  distortion  errors  can  have  several  sources.  They  include: 

•  Aberrations  in  the  electron  optics,  particularly  for  large  scan  fields. 

•  Drift  of  the  electron-optical  axis  due  to  contamination  charging  and  stray  fields. 

•  Drift  in  the  deflection  electronics. 

•  Deviations  from  flatness  and  orthogonality  of  the  interferometer  mirrors  on  the 
workpiece  stage. 

•  Z-axis  deviations  due  to  workpiece  positioning  and  warpage  or  stage  misalign¬ 
ment. 

Some  of  the  distortions  are  reproducible  and  could,  in  principle,  be  compensated  in 
the  software;  others  are  not.  Typical  positioning  specifications  of  e-beam  systems  are 
±  1/10  micron  to  ±  1/4  micron.  There  has  been  little  independent  verification  of 
these  specificatons.  Distortion  errors  are  most  apparent  when  one  attempts  to  draw 
a  pattern  which  extends  over  more  than  one  scan  field.  Lines  which  continue  from 
one  scan  field  to  an  adjacent  scan  field  may  fail  to  butt  smoothly  at  the  boundary  be¬ 
tween  scans  fields,  leading  to  a  pattern  mismatch  referred  to  as  stitching  error. 
When  we  began  our  work,  a  typical  stitching  error  on  the  Honeywell  EBES  was  ±  1 
micron. 

Our  method  of  measuring  e-beam  pattern  distortion  is  the  interferometric  technique 
introduced  by  Wyant  [2-1].  The  principle  of  this  technique  is  simple.  If  an  e-beam  de¬ 
vice  is  used  to  write  a  linear  grating  of  straight,  parallel  lines,  the  resulting  pattern 
can  be  thought  of  as  a  hologram  produced  by  interfering  two  plane  wavefronts.  But, 
because  of  pattern  distortion,  the  lines  are  never  perfectly  straight  nor  equally 
spaced.  A  pattern  of  distorted  straight  lines  can  be  thought  of  as  a  hologram  made 
by  interfering  a  tilted  plane  wave  with  an  aberrated  plane  wave.  Experimentally,  a 
Mach-Zehnder  interferometer  is  set  up  to  evaluate  the  aberrated  plane  wave.  Once 
the  interferogram  is  digitized  and  analyzed,  a  contour  map  for  the  entire  pattern 
showing  lines  of  constant  distortion  in  a  particular  direction  can  be  obtained.  By 
combining  two  contour  maps  for  distortion  in  orthogonal  directions,  the  distortion  in 
any  arbitrary  direction  can  be  found. 


E47069 


2-2 


To  perform  the  e-beam  distortion  analysis,  we  fabricated  linear  gratings  and  rec¬ 
tangular,  crossed  gratings.  A  crossed  grating,  as  shown  in  Figure  2-2,  allows  or¬ 
thogonal  components  of  distortion  to  be  tested  with  a  single  pattern.  When  the 
crossed  grating  is  illuminated  with  a  plane  wave,  several  diffracted  orders  are 
produced  as  shown  in  Figure  2-3.  An  interferometer  arrangement,  as  depicted  in 
Figure  2-4,  was  used  for  interfering  the  +N  and  -N  orders.  Typical  interferograms 
appear  in  Figure  2-5.  A  deviation  from  linearity  of  one  fringe  in  the  interferogram 
represents  a  pattern  distortion  of  1/2N  grating  periods. 

As  indicated  earlier,  our  initial  tests  of  the  Honeywell  EBES  showed  distortion  er¬ 
rors  of  ±  1  micron.  We  took  several  steps  toward  reducing  these  errors: 

•  Scan  field  —  the  size  was  reduced  to  1  mm  x  1  mm  from  1.6  mm  x  1.6mm  to 
minimize  aberrations. 

•  Software  compensation  —  the  functional  form  of  the  distortions  was  fit  to  a 
nine-point  polynomial  and  the  coefficients  were  updated  periodically  during 
drawing  to  allow  for  temporal  variations  in  the  distortion.  The  coefficients  were 
used  to  register  the  pattern  to  substrate  marks  as  shown  in  Figure  2-6. 

•  Hardware  improvements  —  some  sources  of  magnetic  interference  were 
eliminated  and  special  magnetic  shielding  was  installed  to  reduce  the  effects  of 
others. 

These  efforts  brought  us  to  the  point  where  we  could  draw  patterns  with  distortion 
errors  of  ±  0.15  micron  in  the  X-scan  and  ±  0.30  micron  in  the  Y-scan  directions. 


Figure  2-2.  Rectangular  Crossed  Grating  Used  to  Measure  E-Beam  Pattern 
Distortion 
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Figure  2*3.  Diffraction  Pattern  of  the  Crossed  Grating 
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Figure  2-4.  Interferometer  Arrangement  for  Measuring  E-Beam  Pattern 
Distortion  Using  the  Plane  Wave  Interference  Technique 
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(a)  Single  scan  field 


(b)  2X2  multi-scan  field 


Figure  2-5.  Distortion  Tests  of  the  Honeywell  EBES.  These  typical 

interferograms  are  of  linear  gratings  with  a  5-micron  period  and 
1  mm  x  1  mm  scan  field  size. 


Figure  2-6.  Software  Distortion  Correction.  Aberrations  in  the  electron  optics 
are  compensated  by  stretching  the  e-beam  field. 
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Comparison  with  DICOMED  Image  Recorder 


Conventionally,  computer  generated  holograms  have  been  drawn  at  large  scale  on  a 
graphic  plotting  device  and  then  photo-reduced  to  final  size  on  high-resolution  holo¬ 
graphic  film.  Experience  has  shown  that  plotter  distortion  is  one  of  the  larger  error 
sources  in  holograms  produced  by  this  method.  In  this  research  program,  we  sought 
to  measure  the  limitations  of  the  optical  recording  and  photo-reduction  method  of 
hologram  fabrication  in  order  to  assess  the  benefits  of  direct  e-beam  fabrication.  We 
chose  the  DICOMED  Image  Recorder  [2-2]  as  our  optical  recorder  for  this  study. 

Three  crossed  gratings  of  different  spatial  frequencies  were  plotted  by  the 
DICOMED  Image  Recorder  each  month  for  four  consecutive  months  and  then  photo- 
reduced  to  a  final  size  of  12  mm  diameter.  A  distortion  analysis  was  performed  on 
these  crossed  gratings  using  the  same  testing  technique  as  for  e-beam  writing.  Fig¬ 
ure  2-7  shows  a  typical  interferogram  from  the  crossed  grating  test  using  fourth  or¬ 
ders  (N  =  ±  4).  This  crossed  grating  had  200  lines  per  radius  in  each  direction. 


Figure  2-7.  Distortion  Test  of  the  DICOMED  Image  Recorder.  The 

interferograms  were  created  using  the  N  =  ±  4  orders  of  a  crossed 
grating  with  200  lines  per  radius.  They  show  the  components  of 
distortion  in  orthogonal  directions. 


In  this  distortion  analysis,  we  looked  at  three  different  aspects  of  the  distortion  er¬ 
rors: 


•  Distortion  versus  time  (repeatability) 

•  Distortion  versus  polarity  of  the  image 

•  Distortion  versus  photo-reduction  process 
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In  general,  distortions  in  orthogonal  directions  had  similar  values.  Plotting  distor¬ 
tion  did  fluctuate  from  month  to  month  within  the  six-month  time  frame,  but  stayed 
within  the  instrument  specification  of  0.25  percent.  The  dominant  distortion  was 
found  to  be  astigmatism.  Positive  images  had  the  least  error.  Measured  errors  in¬ 
cluded  plotter  distortion  as  well  as  lens  distortion  in  the  photo-reduction  process.  It 
was  found  that  photo-reduction  could  contribute  as  much  as  0.02-percent  distortion 
error  to  the  final  hologram.  The  smaller  the  final  size  of  the  image, ‘the  larger  the  to¬ 
tal  distortion  error. 

We  have  concluded  that  the  net  pattern  distortion  for  the  DICOMED  Image  Re¬ 
corder  is  less  than  0.10  percent,  with  a  usable  resolution  of  2000  elements  per 
diameter.  This  can  be  used  to  encode  wavefronts  with  a  maximum  slope  of  about  125 
waves  per  radius,  for  which  the  maximum  wavefront  error  would  be  0.25  wave.  By 
comparison,  a  1  cm  diameter  e-beam  generated  pattern  with  ±  0.25-micron  accuracy 
has  2  x  104  resolution  elements  per  diameter  and  0.01-percent  distortion.  This  result 
confirms  the  benefits  of  e-beam  fabrication  of  computer  generated  holograms. 

Hologram  Encoding  Method 

For  simple  wavefront  holograms,  only  the  phase  need  be  encoded  since  the 
amplitude  is  constant  over  the  entire  hologram  plane.  Wyant  and  Bennett  [2-3]  plot 
wavefront  holograms  by  drawing  the  centers  of  fringes  formed  by  the  conceptual  in¬ 
terference  between  the  desired  wavefront  and  a  tilted  plane  wave.  In  contrast  to  a 
Fourier  transform  hologram,  the  wavefront  hologram  generally  has  a  simple  func¬ 
tional  form  that  can  be  quickly  evaluated  at  any  point  on  the  hologram  plane. 
Loomis  [2-4]  has  extended  this  fringe-drawing  method  to  calculate  the  phase  of  the 
wavefront  at  every  pixel.  For  patterns  of  up  to  4096  x  4096  pixels,  this  method  uses 
reasonable  amounts  of  computer  time.  For  e-beam  computer  generated  holograms, 
however,  calculation  of  the  wavefront  at  each  pixel  becomes  prohibitive.  We  have 
therefore  modified  Loomis’  method  in  order  to  develop  software  for  producing  prac¬ 
tical  e-beam  computer  generated  wavefront  holograms. 

Our  initial  encoding  method  was  geared  toward  the  pattern  generation  capability  of 
the  Honeywell  EBES.  Since  the  Honeywell  pattern  generator  was  designed  to  effi¬ 
ciently  expose  arbitrarily  oriented  rectangular  areas  of  a  workpiece,  we  chose  to  rep¬ 
resent  a  fringe  by  a  series  of  connected  rectangles  of  various  sizes  and  orientations 
as  depicted  in  Figure  2-8.  A  typical  rectangle  is  perhaps  5  microns  wide  and  25  mi¬ 
crons  long.  Since  our  pixel  size  is  0.3  micron,  this  means  that  one  rectangle  contains 
about  1400  pixels.  Because  only  half  of  the  pixels  in  a  hologram  are  exposed,  we  can 
specify  a  hologram  with  109  pixels  by  giving  the  locations  of  only  3.6  x  105  exposure 
rectangles.  Each  rectangle  is  specified  by  six  numbers  of  16  bits  each.  Therefore,  a 
hologram  can  be  specified  by  a  total  of  3.4  x  107  bits.  Specification  by  Loomis’  meth¬ 
od  would  require  109  bits.  Our  encoding  method  therefore  reduces  the  number  of  bits 
needed  to  specify  a  computer  generated  hologram  by  a  factor  of  about  30. 
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SINGLE  SCAN  FIELD 


Figure  2-8.  Hologram  Encoding  Method  for  Honeywell  EBES.  The  pattern  is 

composed  of  rectangles  in  various  orientations.  Each  rectangle  is  an 
array  of  pixels  on  0.3-micron  centers.  The  exposure  rectangle  is  the 
basic  unit  of  data. 
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The  bulk  of  our  hologram  encoding  is  performed  during  a  preprocessing  step  on 
Honeywell’s  Multics  computer  system.  The  programs  are  written  in  PL1.  A  holo¬ 
gram  is  specified  by  its  wavefront  phase  function  <£(x,y).  A  subroutine  must  exist  for 
evaluating  this  function.  We  are  presently  set  up  to  handle  functions  of  the  form 
[2-5] 


<t>  (x,y)  =  tilt  x  y  +  - •  —  +  a2r2  +  a4r4  +  a6r6  +  a8r8  +  a10r10  +  a12r12 

1+  \1  (1  -  pC^r2) 

where  r2  =  x2  +  (y  -  yj^.  Other  functional  forms  can  easily  be  accommodated  by  re¬ 
placing  the  subroutine. 

Given  a  phase  function  <Mx,y),  the  program  locates  the  centers  of  fringes  as  defined 
by  the  equation  <Mx,v)  =  n,  where  n  assumes  integer  values.  Fringes  are  traced 
across  the  hologram  starting  with  the  fringe  of  lowest  n.  The  centerline  of  each 
fringe  is  approximated  by  a  series  of  connected  straight-line  segments  along  which 
the  phase  4(x,y)  does  not  depart  from  its  integer  value  by  more  than  a  specified 
amount,  typically  0.05.  Since  the  line  segments  are  chosen  to  be  as  long  as  possible, 
this  results  in  longer  segments  where  the  fringes  are  straight,  and  shorter  segments 
where  the  fringes  have  greater  curvature.  The  program  widens  each  line  segment  to 
form  a  rectangle,  chosen  so  as  to  give  a  phase  difference  of  0.5  across  the  width. 

Each  rectangle  is  drawn  as  a  square  array  of  pixels,  usually  on  0.3  micron  centers. 
This  array  is  referred  to  as  an  exposure  rectangle  and  is  the  basic  unit  of  data  ac¬ 
cepted  by  the  e-beam  pattern  generator.  Exposure  rectangles  larger  than  a  max¬ 
imum  size  are  partitioned  into  smaller  rectangles  prior  to  sorting  into  scan  fields.  A 
scan  field  consists  of  all  rectangles  that  are  to  be  drawn  before  repositioning  the 
workpiece.  Rectangles  falling  on  a  boundary  are  arbitrarily  assigned  to  one  of  the 
scan  fields  they  occupy.  After  sorting,  the  exposure  rectangles  are  written  to  mag¬ 
netic  tape  for  reading  by  the  e-beam  computer. 

A  typical  hologram  of  10  mm  diameter  requires  88  scan  fields  and  may  contain  near¬ 
ly  2  x  10s  exposure  rectangles.  It  can  be  drawn  in  about  90  minutes. 

CAMBRIDGE  ELECTRON  BEAM  MICROFABRICATOR  (EBMF) 

In  August  of  1980,  an  EBMF-2  e-beam  system  from  Cambridge  Instruments  was  de¬ 
livered  to  Honeywell’s  Solid  State  Electronics  Division.  Figure  2-9  is  a  photograph  of 
this  system.  We  interferometrically  tested  patterns  of  5-micron  linear  gratings  writ¬ 
ten  by  the  EBMF.  Multi-scan-field  patterns  having  three  different  scan  field  sizes,  1 
mm  x  1  mm,  2  mm  x  2  mm  and  4  mm  x  4  mm,  were  evaluated.  All  patterns  had  the 
same  overall  size  of  8  mm  x  8  mm.  Obvious  distortion,  predominately  pincushion  ef¬ 
fect,  was  observed  in  the  4  mm  x  4  mm  scan  field  pattern  (Figure  2-10).  For  patterns 
drawn  with  a  1  mm  x  1  mm  scan  field  size  (Figure  2-11),  the  pattern  distortion  was 
less  than  ±  0.2  micron.  This  corresponds  to  a  writing  accuracy  of  0.004  percent  over 
a  1-cm  diameter  image. 
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2-9.  Cambridge  EBMF  E-Beam  System.  At  the  far  end  are  the  electron 
optical  column  and  operator's  console.  In  the  foreground  are  the 
system  PDP1 1  computer  and  peripherals. 


Figure  2-10.  Distortion  Test  of  the  Cambridge  EBMF  E-Beam  System  with  4mm 
x  4mm  Scan  Field.  The  8mm  x  8mm  pattern  shows  obvious 
distortion. 


Figure  2-11.  Distortion  Test  of  the  Cambridge  EBMF  E-Beam  System  with  1mm 
x  1mm  Scan  Field.  The  8mm  x  8mm  pattern  shows  only  0.34- 
micron  distortion. 


Because  of  its  superior  writing  accuracy  (even  without  software  compensation)  and 
higher  resolution,  the  Cambridge  EMBF-2  soon  supplanted  the  Honeywell  EBES  for 
writing  of  computer  generated  holograms.  The  Honeywell  system  has  since  been 
abandoned  and  the  Cambridge  system  upgraded  to  an  EBMF-6.  The  remainder  of 
this  section  describes  hologram  fabrication  with  the  EBMF-6. 


Fabrication  Procedure 

Our  current  procedure  for  fabricating  computer  generated  holograms  with  the 
EBMF  involves  several  discrete  steps: 


•  The  hologram  is  encoded  and  written  to  magnetic  tape  in  Cambridge  source  pat¬ 
tern  data  (SPD)  format.  This  computation  is  performed  on  the  Honeywell  Mul- 
tics  computer  system. 

•  The  data  is  “fractured”  into  scan  fields  and  converted  to  binary  pattern  data 
(BPD)  for  drawing.  This  step  utilizes  Cambridge  software  and  runs  on  a  PDP11 
similar  to  the  e-beam  computer. 

•  The  Cambridge  EBMF  is  used  to  expose  the  hologram  pattern  in  e-beam  resist. 

•  After  exposure  and  development,  the  pattern  is  either  wet  etched  to  produce  a 
binary  chrome-on-glass  hologram  or  else  ion  milled  to  produce  an  all-glass 
binary  phase  hologram. 

Hologram  Encoding  Methods 

Our  switch  to  the  Cambridge  EBMF  for  hologram  writing  required  the  development 
of  new  encoding  algorithms.  The  reason  for  this  is  that,  while  the  Honeywell  EBES 
could  draw  rectangles  of  any  orientation,  the  Cambridge  machine  is  only  capable  of 
drawing  rectangles  at  a  few  discrete  angles.  Consequently,  we  have  developed  sever¬ 
al  algorithms  for  reducing  a  hologram  to  a  set  of  rectangles  with  horizontal  and  ver¬ 
tical  sides.  These  algorithms,  described  below,  can  be  compared  on  a  basis  of  number 
of  wavefront  evaluations  and  amount  of  data  storage  required  per  hologram  pixel. 
Here  the  term  pixel  refers  to  resolution  actually  utilized,  which  can  be  considerably 
less  than  the  full  distortion-free  resolution  available.  Frequently,  pixel  size  will  dif¬ 
fer  in  the  two  axes. 

Poly  Format  —  Our  first  encoding  algorithm  for  the  Cambridge  EBMF  is  a  direct 
adaptation  of  the  previously  described  fringe-drawing  method.  Fringes  are  approx¬ 
imated  by  a  series  of  points  connected  by  straight-line  segments.  Standard  Cam¬ 
bridge  data-processing  programs  are  then  used  to  fracture  the  polygonal  fringes  into 
primitive  shapes  [2-6].  Poly  format  results  in  efficient,  compact  source  code,  as  did 
the  earlier  algorithm  from  which  it  derives.  Typical  performance  is  0.008  wavefront 
evaluations  and  0.4  bits  of  source  pattern  data  per  hologram  pixel.  However,  for  any 
but  the  simplest  holograms,  subsequent  fracturing  consumes  enormous  amounts  of 
computer  time  and  produces  prodigious  quantities  of  binary  pattern  data  that  quick¬ 
ly  fill  all  available  storage  devices. 

Staircase  Format  —  This  encoding  method  can  best  be  described  by  referring  to 
Figure  2-12.  Fringes  are  approximated  by  a  number  of  rectangles  with  horizontal 
and  vertical  sides.  The  resulting  “staircase”  boundary  is  quite  visible  in  regions 
where  the  fringes  lie  at  significant  angles  to  the  scan  axes.  The  staircase  can  be 
made  as  fine  as  0.1  micron  if  one  is  willing  to  accept  the  penalty  in  size  of  data  set. 
However,  it  is  sufficient  for  the  step  size  to  be  a  small  fraction  (typically  0.05)  of  the 
fringe  period,  and  steps  are  therefore  made  as  large  as  possible  while  still  satisfying 
this  criterion.  It  is  this  ability  to  use  large  steps,  and  hence  fewer  rectangles,  which 
constitutes  the  principle  advantage  of  Staircase  format.  Fracturing  to  binary  pattern 
data  is  quicker  than  for  Poly  format  and  actually  reduces  the  size  of  data  set.  Defin¬ 
ing  pixel  size  as  the  scale  of  the  staircase  boundary  (variable  across  the  hologram), 
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Figure  2-12.  STAIRCASE  Encoding  Method  for  Cambridge  EBMF.  The  rectangle 
is  the  basic  unit  of  data.  Only  four  rectangle  orientations  are 
’  allowed,  hence  the  stepped  pattern. 

t 

•  Staircase  format  requires  about  0.1  wavefront  evaluation  and  12  bits  of  source  pat- 

I  tern  data  per  hologram  pixel. 

‘  Raster  Format  —  The  motivation  for  this  third  encoding  method  is  to  produce  frac- 

!  tured  binary  pattern  data  directly  and/or  to  produce  data  for  our  raster-scan  ma- 

*  chine.  The  Raster  encoding  method  abandons  fringe  tracing  and  instead  proceeds 

!  scan  field  by  scan  field  and  pixel  by  pixel.  Pixel  size  is  fixed  on  the  basis  of  space- 

1  bandwidth  requirements,  but  generally  exceeds  the  0.1  micron  e-beam  address  reso¬ 

lution.  A  typical  pixel  size  is  0.4  micron  by  3.2  microns.  The  Raster  algorithm  is  not 
yet  fully  implemented,  in  that  wavefront  evaluation  occurs  at  every  pixel  and  the 
output  is  still  in  Cambridge  SPD  format.  Because  of  its  fixed  pixel  size,  this 
algorithm,  even  when  using  run-length  encoding,  produces  somewhat  more  source 

j  pattern  data  than  the  Staircase  algorithm. 

Present  Capabilities  and  Limitations 

The  capabilities  and  limitations  of  the  Cambridge  EBMF  as  a  hologram  writing  tool 
,  are  summarized  below.  Staircase  encoding  format  is  assumed. 

I 
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Linewidth  and  Pattern  Accuracy  —  At  present,  we  are  capable  of  routinely  fabri¬ 
cating  e-beam  computer  generated  holograms  with  spatial  accuracy  (distortion-free 
resolution)  of  0.4  microns.  Line  widths  of  1  micron  are  readily  achievable,  and  nar¬ 
rower  lines  are  probably  not  useful  without  improved  pattern  accuracy. 

Wavefront  Specification  —  The  Staircase  encoding  algorithm  will  accept  any 
wavefront  $(x, y)  having  sufficient  tilt  to  assure  d^x.yVdy  >  0. 

Maximum  Hologram  Size  —  Hologram  size  is  essentially  data-limited.  Staircase 
format  allows  about  4  x  105  exposure  rectangles  per  reel  of  magnetic  tape.  This  is 
enough  to  encode  holograms  having  4  x  lCr/A <t>  pixels,  where  &<t>  is  the  desired 
wavefront  accuracy  (A<£  =  0.05  typical).  A  typical  hologram  is  10  to  12  mm  in 
diameter,  contains  300  to  400  fringes,  and  fills  most  of  1  reel  of  tape.  Workstage 
travel  places  an  absolute  upper  limit  of  100  mm  x  100  mm  on  hologram  size. 

Diffraction  Efficiency  —  Efficiencies  near  the  theoretical  maximum  of  10.1  per¬ 
cent  for  binary  absorption  holograms  and  40.5  percent  for  binary  phase  holograms 
are  attainable. 

Fabrication  Time  —  Turnaround  time  is  limited  primarily  by  e-beam  availability 
and  varies  from  1  day  to  several  weeks.  Actual  fabrication  time  is  typically  30 
minutes  to  encode,  several  hours  to  fracture,  45  minutes  to  expose,  and  30  minutes 
to  develop  and  ion  mill.  The  exposure-limited  writing  speed  of  the  EBMF  is  0.02  to 
1.0  cm2/min,  depending  on  required  resolution. 

Potential  for  Further  Improvement 

Significant  potential  exists  for  increasing  the  hologram  writing  capability  of  the 
EBMF  e-beam  through  software  improvements.  At  present  we  are  severely  data-lim¬ 
ited.  Cambridge  SPD  format  requires  about  12  bits  per  hologram  pixel;  this  is  clear¬ 
ly  excessive.  Fracturing  time  could  be  eliminated  and  data  quantity  halved  by  gen¬ 
erating  BPD  format  data  directly  on  Multics.  Generating  data  directly  on  the  e- 
beam's  PDPll  as  exposure  proceeds  might  eliminate  the  data  bottleneck  entirely. 
Software  distortion  correction  could  improve  the  spatial  accuracy  to  0.2  microns. 

The  Cambridge  e-beam’s  pattern  generation  capability  is  not  well  suited  to  the 
fringe  drawing  method  of  hologram  generation.  In  general,  commercially  available 
vector-scan  e-beam  systems  are  aimed  toward  sparse,  rectilinear  patterns,  while 
raster-scan  systems  sacrifice  resolution  (0.5  to  1.0  micron  typical)  for  drawing  speed. 

If  the  data  problem  was  overcome,  perhaps  by  converting  to  raster  scanning  and 
generation  of  data  while  drawing,  a  higher  resolution  e-beam  such  as  the  EBMF 
could  likely  provide  0.2-micron  spatial  accuracy  over  a  10  cm  x  10  cm  area  with  a 
drawing  speed  of  1  cm2/min.  This  would  correspond  to  a  space-bandwidth  product,  or 
distortion-free  pixel  count,  of  2.5  x  1011  which  compares  favorably  with  in- 
terferometrically  recorded  holograms. 
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Section  3 
Applications 


Our  second  main  objective,  after  having  demonstrated  the  feasibility  of  e-beam  com¬ 
puter  generated  holograms,  has  been  to  demonstrate  useful  applications  of  these  op¬ 
tical  elements.  We  have  demonstrated  applications  in  two  major  areas:  optical  tes¬ 
ting  and  optical  computing. 

We  have  developed  an  aspheric  optical  testing  capability  utilizing  e-beam  computer 
generated  holograms.  This  technology  has  been  used  to  provide  routine  testing  of 
aspheric  diamond-turned  optics  at  our  Electro-Optics  Division.  We  have  also  de¬ 
veloped  a  scheme  for  optical  vector-matrix  multiplication  which  utilizes  an  e-beam 
generated  holographic  mask.  This  approach  results  in  a  fairly  simple  optical  system 
offering  both  large  numerical  range  and  high  accuracy.  Simple  holographic  masks 
have  been  fabricated  and  tested. 

In  the  remainder  of  this  section,  we  describe  in  detail  our  application  of  e-beam  com¬ 
puter  generated  holograms  to  both  aspheric  optical  testing  and  optical  vector-matrix 
multiplication. 

ASPHERIC  OPTICAL  TESTING 

The  high  performance  and  low  cost  requirements  of  modern  optical  systems  demand 
the  frequent  use  of  aspheric  surfaces  and  thus  the  ability  to  test  these  surfaces. 
Often,  aspheric  testing  is  impossible  when  conventional  master  surfaces  are  used  as 
test  plates.  Therefore,  deep  aspheric  surfaces  are  usually  avoided  by  optical  de¬ 
signers.  This  problem  can  be  alleviated  by  storing  the  desired  aspheric  wavefront  in 
the  form  of  a  computer  generated  hologram  [3-1,  3-2,  3-3].  This  stored  wavefront  can 
later  be  used  in  the  interferometric  testing  of  supposedly  identical  aspheric  optical 
elements,  much  as  a  test  plate  is  used  in  the  conventional  testing  of  spherical  optical 
elements. 

The  Honeywell/Tropel  interferometer  used  in  this  work  is  a  custom  designed  instru¬ 
ment  capable  of  testing  aspheric  surfaces  using  computer  generated  holograms.  Fig¬ 
ure  3-1  shows  the  physical  appearance  of  this  interferometer.  Figure  3-2  is  an  ex¬ 
ploded  view  depicting  all  the  optical  elements.  Basically,  the  instrument  is  a 
Twyman-Green  interferometer  in  which  both  the  aspheric  test  beam  and  the 
spherical  reference  beam  pass  through  the  hologram.  The  first-order  reference  beam 
and  zero-order  (undiffracted)  test  beam  are  then  combined  to  produce  an  interference 
pattern.  Much  of  the  apparent  complexity  results  from  the  use  of  high-efficiency 
polarizing  beam  splitters.  A  detailed  description  of  this  interferometer  is  given 
elsewhere  [3-4]. 

To  illustrate  the  procedure  for  using  e-beam  computer  generated  holograms  to  test 
aspheric  optics,  we  use,  as  an  example,  the  testing  of  a  4-inch  diameter  f/2  parabolic 
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Figure  3-1 .  A  Photograph  of  the  Honey  well/Tropel  Holographic  and  Shearing 
Interferometer  in  Operation 

mirror  [3-5].  A  second  example  of  a  tilted  plane-parallel  plate  will  then  demonstrate 
the  capability  of  testing  deep  aspheric  optics  lacking  rotational  symmetry. 

Example  1 :  4-inch  Diameter  f/2  Parabola 

The  procedures  involved  in  the  testing  of  the  f/2  parabola  are  as  follows: 

Step  1  —  A  mathematical  description  of  the  f/2  parabola  is  used  in  a  ray  tracing  pro¬ 
gram  along  with  the  optical  design  data  for  the  interferometer.  The  computer  calcu¬ 
lates  the  best  position  for  the  parabola,  and  also  the  fringe  pattern  that  would  be 
produced  at  the  hologram  plane  if  a  perfect  test  wavefront  were  interfered  with  a 
reference  wavefront  at  a  particular  non-zero  incident  angle.  Known  imperfections  in 
the  interferometer  optics  can  be  compensated  for  in  the  ray  tracing. 

Step  2  —  The  mathematical  description  of  the  fringe  pattern  is  used  to  encode  and 
fabricate  an  e-beam  computer  generated  hologram  using  the  methods  discussed  in 
the  previous  section.  The  finished  hologram  and  its  diffraction  pattern  are  shown  in 
Figures  3-3  and  34. 
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Figure  3-4.  Diffraction  Pattern  Due  to  the  Hologram  of  Figure  3-3. 


E47069 


3-5 


Step  3  —  When  the  finished  e-beam  computer  generated  hologram  is  inserted  prop¬ 
erly  at  the  hologram  plane  of  the  interferometer  system,  both  the  test  beam  and  the 
reference  beam  are  diffracted  by  the  hologram  into  plus  and  minus  first-order 
beams.  When  the  zero-order  test  beam  and  the  zero-order  reference  beam  are  super¬ 
imposed,  the  02  parabola  is  tested  without  aid  of  the  computer  generated  hologram, 
and  the  interferogram  of  Figure  3-5a  results.  It  shows  about  40  waves  of  spherical 
aberration.  This  Twyman-Green  interferogram  would  be  exceedingly  difficult  to 
analyze  in  detail.  However,  its  complexity  is  incorporated  into  the  hologram,  as 
shown  by  the  interferogram  in  Figure  3-5b.  This  latter  interferogram,  produced  by 
interfering  zero  and  first  orders  of  a  reference  beam,  simulates  the  Twyman-Green 
test  of  a  perfect  f/2  parabola  with  compensation  for  known  instrumental  errors. 
When  the  reference  beam  is  properly  adjusted,  the  first-order  reference  and  the  zero- 
order  test  beams  can  be  superimposed  to  produce  the  interferogram  of  Figure  3-5c. 
This  interferogram  gives  the  residual  wavefront  error  of  the  f/2  parabola.  If  the  pa¬ 
rabola  were  perfect,  the  fringes  would  be  straight,  parallel  and  equally  spaced.  For 
comparison,  the  results  of  an  autocollimation  test  are  shown  in  Figure  3-6.  The  two 
tests  are  in  good  agreement. 

Example  2:  Asymmetric  Deep  Aspheric 

The  f/2  parabola  test  demonstrated  the  feasibility  of  using  e-beam  computer  gener¬ 
ated  holograms  to  test  known  conic  surfaces  of  moderate  asphericity.  For  our  second 
example,  we  designed  an  asymmetric  deep  aspheric  wavefront  using  a  concave  mir¬ 
ror  and  a  tilted  plane-parallel  plate  as  the  test  surface.  Figure  3-7  shows  the  ex¬ 
perimental  setup  used  to  obtain  this  wavefront.  The  amount  of  aspheric  aberration 
in  this  wavefront  can  be  selected  by  tilting  the  plane-parallel  plate  placed  between 
the  concave  mirror  the  the  diverger.  Table  3-1  lists  wavefront  aberration  as  a  func¬ 
tion  of  parallel  plate  tilt  angle. 

Based  on  the  coefficients  obtained  from  ray  tracing,  we  fabricated  e-beam  computer 
generated  holograms  to  test  this  aspheric  wavefront.  This  required  a  departure  from 
our  normal  parametric  description,  but  the  change  was  easily  accommodated  be¬ 
cause  of  the  modularity  of  our  encoding  software.  Holograms  were  generated  cor¬ 
responding  to  parallel  plate  tilt  angles  of  20  and  25  degrees.  The  rather  unusual  dif¬ 
fraction  patterns  of  these  holograms  are  shown  in  Figures  3-8  and  3-9. 

Figure  3-10  is  a  Twyman-Green  interferogram  of  the  test  wavefront  with  a  plate  tilt 
angle  of  20  degrees  (analogous  to  Figure  3-5a  for  the  f/2  parabola).  Figure  3-11  shows 
the  e-beam  computer  generated  hologram  test  of  this  asymmetric,  aspheric 
wavefront.  This  test  clearly  demonstrates  the  capability  of  testing  generalized 
aspheric  surfaces  with  e-beam  computer  generated  holograms. 
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Figure  3-5.  Testing  of  the  f/2  Parabola  with  an  E-Beam  Computer  Generated 
Hologram 


Figure  3-6.  Autocollimation  Test  of  the  f/2  Parabola.  Surface  deviation  at  6328$, 
is  0.035  waves  rms,  0.171  waves  pv. 
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Figure  3-7.  Experimental  Arrangement  to  Obtain  Asymmetric  Aspheric 
Wavefront 
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Table  3-1.  Asymmetric  Aspheric  Wavefronts 


I 


Tilt  Angle 

Slope  (Waves/Radius) 

15* 

71 

17.3° 

83 

20* 

97 

25° 

123 

Figure  3-8.  Diffraction  Pattern  of  E-Beam  Computer  Generated  Hologram  to 

Simulate  Asymmetric  Aspheric  Wavefront  for  Plate  Tilt  Angle  of  20 
Degrees 


Figure  3-9.  Diffraction  Pattern  of  an  E-Beam  Computer  Generated  Hologram  to 
Simulate  Asymmetric  Aspheric  Wavefront  for  Plate  Tilt  Angie  of  25 
Degrees 
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Figure  3-10.  Twyman-Green  Interferogram  of  Asymmetric  Aspheric  Wavefront 
for  Plate  Tilt  Angle  of  20  Degrees 


Figure  3-11.  E-Beam  Computer  Generated  Hologram  Test  of  Asymmetric 

Aspheric  Wavefront  for  Plate  Hit  Angle  of  20  Degrees.  Surface 
deviation  at  6328&  is  0.04  waves  rms,  0.29  waves  pv. 


OPTICAL  VECTOR-MATRIX  MULTIPLICATION 

There  has  recently  been  considerable  interest  in  optical  computing  since  it  offers 
very  high  computation  throughput  rates  for  mathematical  operations  amenable  to 
parallel  computation.  One  class  of  such  operations,  vector-matrix  multiplication 
with  a  fixed  matrix,  can  be  used  for  performing  discrete  Fourier  transforms,  coordi¬ 
nate  transformations,  pattern  classification,  and  many  other  computations  [3-6,  3-7, 
3-8, 3-9, 3-10].  The  general  matrix-vector  multiplication  may  be  written  as: 

N 

ym  "  2  Hmn  x“  (m’1’  2>  •  •  -W 

n-l 
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Our  optical  approach  to  performing  this  computation  uses  N  light  sources  to  repre¬ 
sent  the  components  of  xn  of  the  input  vector,  M  detectors  to  represent  the  compo¬ 
nents  of  ym  of  the  output  vector,  and  suitable  optics  to  assure  that  a  fraction  Hmn  of 
the  light  from  source  xn  gets  to  each  detector  ym.  The  problem  can  be  suitably  scaled 
so  that  all  parameters  fall  within  acceptable  positive  ranges.  Optics  to  perform  the 
function  of  the  matrix  H  remain  fixed,  while  the  N  sources  are  modulated  to  repre¬ 
sent  various  input  vectors  x. 

In  principle,  the  performance  of  this  optical  computer  is  dependent  on  a  number  of 
considerations  involving  the  optics,  detectors  and  sources.  In  practice,  accuracy  is 
often  limited  by  matrix  mask  imperfections,  while  speed  is  limited  by  the  amount  of 
light  reaching  the  detectors.  For  this  reason,  our  work  has  focused  on  efficient  optics 
to  precisely  distribute  light  among  the  various  detectors. 

In  most  schemes  for  optical  vector-matrix  multiplication,  the  matrix  is  encoded  as  a 
rectangular  array  of  apertures  or  gray  tones  in  a  mask.  This  approach,  represented 
in  Figure  3-12,  encounters  several  limitations.  A  complicated  optical  system  is  re¬ 
quired  in  order  to  illuminate  and  receive  light  from  specific  columns  and  rows  of  the 
matrix  mask.  Much  light  is  discarded  in  providing  uniform  illumination  to  the 
mask,  with  the  mask  passing  only  about  half  of  what  remains.  Numerical  range  and 
accuracy  are  limited  by  the  space-bandwidth  product  of  the  mask  (generally  less 
than  106  with  conventional  plotting  techniques).  Small  matrix  elements  result  in 
small  apertures  with  low  relative  accuracies.  If  results  differ  from  those  intended,  it 
is  generally  difficult  to  modify  a  mask  except  by  starting  anew. 
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Figure  3-12.  A  general  scheme  for  Optical  Vector-Matrix  Multiplication 


PCGH  Configuration. . 

Ll  Our  approach  to  optical  vector-matrix  multiplication,  depicted  in  Figure  3-13,  is 

based  upon  an  e-beam  generated  diffractive  mask  which  we  will  call  a  partitioned 
computer  generated  hologram  or  PCGH  [3-10].  Each  of  N  PCGHs  is  illuminated  by 
collimated  light  from  a  single  element  of  the  source  array  and  thus  represents  one 
column  of  the  NxM  matrix  mask  in  Figure  3-12.  Each  PCGH  is  partitioned  into  M 
linear  gratings,  which  diffract  light  to  the  M  detectors.  The  optical  power  diffracted 
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Figure  3-13.  Partitioned  Computer  Generated  Hologram  (PCGH),  Vector-Matrix 
Multiplier.  The  several  facets  of  each  PCGH  contain  linear  gratings 
that  deflect  portions  of  the  light  to  the  various  detectors. 

by  a  particular  grating  is  made  proportional  to  the  value  of  the  required  matrix  ele¬ 
ment.  The  PCGHs  are  fabricated  as  binary  chrome-on-glass  holograms  with  the  grat¬ 
ing  pattern  delineated  via  e-beam  lithography. 

Our  e-beam  PCGH  vector-matrix  multiplication  scheme  has  several  advantages  over 
the  scheme  in  Figure  3-12.  E-beam  lithography  offers  a  higher  space-bandwidth 
product,  which  can  translate  into  greater  numerical  accuracy.  Also,  the  PCGHs  are 
in  the  Fourier  plane  of  the  tranform  lens  with  respect  to  the  detectors.  This  means 
that  the  only  requirement  for  light  to  reach  a  particular  detector  is  that  it  be  trav¬ 
eling  in  the  proper  direction  upon  leaving  the  PCGH.  Therefore,  each  input  module, 
consisting  of  source,  collimating  lens  and  PCGH,  may  be  located  anywhere  within 
the  aperture  of  the  transform  lens.  This  same  immunity  to  translation  allows  a 
PCGH  to  be  partitioned  into  facets  in  any  manner  consistent  with  dividing  up  the 
available  light  amongst  the  various  detectors  (providing,  of  course,  that  the  facets  do 
not  become  too  small).  Other  advantages  relate  to  optical  efficiency.  Except  for  dif¬ 
fraction  efficiency  losses,  all  light  striking  the  PCGH  can  be  used.  Light  need  not  be 
wasted  in  achieving  uniform  illumination;  non-uniform  illumination  is  acceptable  so 
long  as  its  effects  are  accounted  for  in  the  partitioning.  Small  facets  associated  with 
lesser  outputs  can  be  made  physically  larger  by  placing  them  where  PCGH  illumina¬ 
tion  is  lowest.  The  various  considerations  which  affect  the  design  of  a  PCGH  are  dis¬ 
cussed  in  the  following  sections. 

PCGH  Diffraction  Analysis 

Figure  3-14  illustrates  a  simple  PCGH  intended  to  produce  10  equal  intensity  out¬ 
puts  when  uniformly  illuminated.  This  PCGH  has  10  equal  area  facets,  each  con¬ 
taining  a  linear  grating  to  diffract  incident  light  to  the  appropriate  detector.  The 
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Figure  3-14.  PCGH  Designed  for  Uniform  Illumination  and  10  Equal  Intensity 
Outputs.  Facets  are  arranged  symmetrically  to  provide  immunity 
to  beam  wander. 

numbers  in  each  facet  indicate  the  channel  assignments.  Spatial  frequencies  of  the 
linear  gratings  are  dictated  by  the  system  geometry.  First-order  diffracted  light  from 
a  facet  of  spatial  frequency  v  will  be  focused  in  the  detector  plane  a  distance  v\  F 
from  the  transform  lens  axis  (Figure  3-13),  where  X  is  the  wavelength  and  F  the 
transform  lens  focal  length.  Our  design  is  for  a  10-element  linear  detector  array. 
This  requires  10  equally  spaced  grating  frequencies.  The  widest  possible  detector 
separation  is  achieved  with  grating  frequencies  (n  +  9)Aj»,  where  Av  is  the  frequency 
separation  and  n  =  1,2,  . .  .,10.  Then  the  unwanted  harmonic  frequencies  from  the 
imperfect  square-wave  gratings  begin  at  20Av  and  do  not  coincide  with  desired  out¬ 
puts. 

Matrix  values  are  encoded  into  the  PCGH  via  grating  area  modulation.  The  holo¬ 
gram  must  therefore  be  divided  into  facets  in  such  a  way  that  the  amount  of  light 
diffracted  by  a  facet  to  its  detector  is  proportional  to  the  required  matrix  element. 
Various  considerations  lead  us  to  partition  the  PCGH  into  facets  along  a  rectangular 
grid. 

Mathematically,  the  transmittance  of  a  facet  can  be  regarded  as  the  product  of  its 
aperture  and  an  infinite  linear  grating.  By  the  Fourier  convolution  theorem,  the  dif¬ 
fraction  pattern  of  this  facet  is  the  diffraction  pattern  of  its  aperture  convolved  with 
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the  delta  function  from  its  grating.  In  other  words,  the  effect  of  the  linear  grating  is 
to  shift  the  location  of  the  diffraction  pattern  produced  by  the  facet  aperture. 

Figure  3-15  indicates  the  diffraction  pattern  due  to  a  square  aperture  of  dimension 
D.  The  main  lobe  has  a  width  of  2XF/D  and  contains  81.5  percent  of  the  energy  pass¬ 
ing  through  the  aperture.  The  sidelobes  form  a  rectangular  array  with  sidelobe 
energy  diminishing  inversely  as  the  square  of  the  distance  from  either  axis  as  in¬ 
dicated  by  Figure  3-16.  Crosstalk  therefore  diminishes  most  rapidly  for  detectors 
spaced  along  a  line  oriented  at  45  degrees  to  the  facet  boundaries.  This  requires  the 
grating  fringes  to  run  diagonally  within  each  facet. 

As  indicated  earlier,  the  separation  of  detectors  is  ApXF.  As  facets  are  made  smaller, 
their  diffraction  patterns  become  larger,  producing  undesirable  crosstalk.  For  this 
reason,  we  impose  a  minimum  square  facet  size  D  to  produce  a  main  diffraction  lobe 
whose  (diagonal)  dimension  (2V2)XF/D  is  2/3  the  (diagonal)  detector  separation  ApXF. 
In  other  words,  the  minimum  facet  size  is  D  =  (3\/2)/Av.  We  size  our  detector 
aperture  to  capture  only  the  main  lobe  of  this  minimum  facet  diffraction  pattern  for 
a  calculated  worst  case  crosstalk  of  37  dB.  A  larger  aperture  would  capture  lesser 
lobes  of  the  channel  of  interest,  but  also  some  greater  lobes  of  adjacent  channels, 
thereby  degrading  the  channel  separation.  Also,  having  zero  intensity  at  the  edges 
of  the  detector  aperture  serves  to  ease  aperture  alignment  tolerances. 

The  above  discussion,  of  course,  considers  crosstalk  when  each  facet  is  a  minimum 
square  facet.  This  would  be  a  severe  constraint  on  system  numerical  range.  In  prac¬ 
tice,  we  form  the  facet  for  each  channel  from  many  of  these  minimum-sized  modular 
subfacets.  Therefore,  the  diffraction  pattern  for  each  channel  is  the  coherent  sum  of 
all  subfacets  for  that  channel.  Crosstalk  can  be  minimized  by  clustering  each  chan¬ 
nel’s  subfacets  into  one  or  two  large  facets,  as  was  done  in  the  PCGH  of  Figure  3-14. 
By  requiring  that  large  facets  be  built  up  of  modular  subfacets,  we  ensure  that,  no 
matter  how  complicated  the  diffraction  pattern  within  a  detector  aperture,  the  in¬ 
tensity  will  still  fall  to  zero  at  the  aperture  edges.  Also,  since  subfacet  boundaries 
are  diagonal  to  the  grating  fringes,  any  bright  sidelobes  will  occur  diagonally  away 
from  adjacent  detectors.  This  phenomenon  is  easily  seen  in  Figure  3-20. 

Partitioning  the  PCGH 

Partitioning  is  the  task  of  assigning  subfacets  to  each  channel  in  such  a  way  as  to 
achieve  the  correct  relative  intensities.  This  task  is  greatly  simplified  by  our  de¬ 
cision  to  adopt  modular  subfacets.  The  energy  incident  on  a  rectangular  facet  is  easi¬ 
ly  computed  for  either  uniform  or  Gaussian  illumination.  Each  subfacet  is  assigned 
entirely  to  one  channel.  We  assume  that  each  channel  receives  several  subfacets 
which  can  be  grouped  into  one  or  more  larger  facets.  Most  of  the  diffracted  energy 
from  one  of  these  larger  facets  should  fall  within  the  detector  aperture.  Therefore, 
the  contribution  of  each  subfacet  to  the  total  detected  energy  is  approximately  the 
product  of  its  incident  illumination  and  its  diffraction  efficiency.  We  need  only  com¬ 
pute  this  energy  for  each  subfacet  and  then  decide  which  subfacets  are  to  be  as¬ 
signed  to  which  channels.  The  error  introduced  by  this  approximation  varies  from 
near  zero  (many  subfacets/facet)  to  18.5  percent  (one  subfacet/facet). 
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Figure  3*16.  Relative  Intensity  of  Sidelobes  (dB)  for  Square  Aperture 

Diffraction.  A  square  detector  captures  the  main  lobe.  Adjacent 
identical  detectors  (dashed  boxes)  must  be  separated  diagonally  to 
avoid  bright  on-axis  sidelobes.  Crosstalk  (dB)  is  indicated  for  three 
alternative  adjacent  detector  locations. 


To  have  the  greatest  flexibility  in  partitioning  the  PCGH,  and  to  minimize  the  ap¬ 
proximation  error,  we  would  like  the  modular  subfacets  to  be  as  small  and  numerous 
as  possible.  Since  they  are  of  size  D  =  (3V2 )/Av,  we  want  a  large  spatial  frequency 
separation  Av.  However,  a  large  spatial  frequency  separation  implies  large  spatial 
frequencies  and  hence  small  grating  periods.  If  the  grating  period  becomes  com¬ 
parable  to  the  e-beam  spot  size,  considerable  grating  duty  cycle  errors,  with  cor¬ 
responding  diffraction  efficiency  errors,  will  result.  For  these  reasons  we  have 
elected  to  use  grating  frequencies  of  60  to  114  lp/mm  (measured  along  either  axis) 
and  0.6  mm  subfacets.  This  gives  us  400  subfacets  in  a  1  cm  x  1  cm  PCGH  for  a 
maximum  numerical  accuracy  of  about  20  dB.  With  Gaussian  illumination,  the  ef¬ 
fective  number  of  subfacets  is  more  than  5000,  since  the  corner  facets  will  have 
about  6  percent  the  intensity  of  illumination  of  the  central  facets,  assuming  a  1  cm 
(1/e2)  beam  diameter. 
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A  well  partitioned  PCGH  can  be  characterized  by  the  following  four  criteria: 

•  Subfacets  for  each  channel  should  add  up  to  the  correct  total  power. 

•  Subfacets  for  each  channel  should  form  a  single,  compact  facet. 

•  In  the  case  of  Gaussian  illumination,  facets  for  the  dimmer  channels  should  be 
located  near  the  edges  so  that  they  may  be  made  larger. 

•  In  the  case  of  Gaussian  illumination,  subfacets  for  each  channel  should  be  locat¬ 
ed  symmetrically  about  the  center  to  provide  some  immunity  to  beam  wander. 

It  is  generally  not  possible  to  simultaneously  satisfy  each  of  these  criteria.  For  ex¬ 
ample,  a  facet  cannot  be  symmetric  about  the  center,  located  near  the  edge,  and  also 
compact.  Some  compromise  is  necessary. 

The  partitioning  algorithm  which  we  have  developed  uses  a  merit  function  to  decide 
if  a  particular  change  to  the  PCGH  is  advantageous.  The  merit  function,  which  em¬ 
bodies  our  four  criteria,  is 

MF  =  CjxE  Irelative  intensity  error| 

+  C2  x  [number  of  facet  corners  +  2  x  facet  perimeter] 

+  C3  x  L  (deviation  from  mean  of  facet  size) 

+  C4  x  [number  of  asymmetrically  assigned  subfacets] 

where  the  summations  are  over  the  various  channels.  The  units  of  facet  perimeter 
and  facet  size  are  D  and  D2  respectively.  The  constants  Cx,  C2,  C3  and  C4  are  ad¬ 
justed  to  give  desired  emphasis  to  the  four  criteria.  This  merit  function  is  a  negative 
trait  in  the  sense  that  we  wish  to  minimize  it. 

While  the  merit  function  provides  a  means  of  distinguishing  the  better  of  two  parti¬ 
tionings  of  a  PCGH,  it  does  not  tell  us  which  partitionings  to  consider.  With  400  sub¬ 
facets  and  10  channels,  there  are  IO400  possible  partitionings!  An  exhaustive  search 
is  out  of  the  question.  We  have  developed  two  algorithms  for  partitioning.  One 
algorithm  takes  a  starting  PCGH  (which  may  be  blank)  and  considers  all  changes  in¬ 
volving  only  one  subfacet.  Changes  which  would  divide  existing  facets  are  rejected. 
Any  other  change  which  reduces  the  merit  function  is  implemented.  Although  this 
algorithm  can  be  somewhat  restrictive,  it  has  tended  to  work  reasonably  well.  A  sec¬ 
ond  algorithm,  which  we  have  also  found  useful,  allows  two  subfacets  to  swap  their 
channel  assignments.  As  is  the  case  with  most  optical  design  codes,  the  operator 
must  guide  the  partitioning  to  some  extent,  redefining  the  merit  function  according 
to  how  the  PCGH  is  developing  and  choosing  which  algorithm  to  apply. 

Figure  3-17  shows  the  partitioning  algorithm  at  work  on  a  PCGH  having  20  x  20  or 
400  subfacets.  The  goal  was  a  hologram  that  would  produce  10  equal  intensity  out¬ 
puts  when  uniformly  illuminated  (40  subfacets/channel).  Facet  boundaries  are  in¬ 
dicated  following  the  first,  one  intermediate  and  the  last  iteration.  Figure  3- 17c 
should  be  compared  to  the  PCGH  of  Figure  3-14  which  included  a  symmetry 
criterion  in  its  merit  function. 
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a)  After  the  first  iteration,  b)  Several  iterations  later,  c>  Final  solution.  The 

facets  branch  from  their  most  facets  are  algorithm  could  find  no 

starting  points  along  the  considerably  more  further  improvements, 

upper  left  margin.  compact. 

Figure  3-17.  Partitioning  of  a  PCGH  for  Uniform  Illumination  and  10  Equal 
Intensity  Outputs 


Trimming  Methods 

To  extend  our  numerical  accuracy  beyond  about  20  dB,  it  is  necessary  to  employ  a 
separate  lithography  step  to  adjust  the  relative  amplitudes  of  the  outputs.  The  abili¬ 
ty  to  do  this  trimming  is  one  of  the  chief  advantages  of  e-beam  lithography  and  the 
PCGH.  There  are  several  possible  methods  of  trimming  a  PCGH  once  it  has  been 
made  and  tested.  The  best  trimming  method  is  to  add  a  negative  facet,  that  is  to  say, 
a  facet  exactly  out  of  phase  with  the  existing  facet.  This  is  easier  than  it  might 
sound:  the  negative  facet  can  be  written  in  space  already  occupied  by  the  existing 
positive  facet,  leaving  a  completely  open  area  in  a  portion  of  the  positive  facet.  Line- 
width  problems  are  overcome  and  registration  requirements  reduced  since,  if  we 
write  over  the  whole  area  to  be  trimmed  (i.e.,  not  just  when  there  is  chrome),  then 
the  phase  and  duty  cycle  of  the  negative  facet  exactly  complement  the  existing 
positive  facet. 

Experimental  Results 

The  experimental  setup  for  demonstration  and  evaluation  of  PCGHs  is  shown  in  Fig¬ 
ure  3-18.  The  setup  resembles  Figure  3-13  except  f'-r  the  detector  array,  which  has 
been  replaced  by  a  single,  movable  detector.  Light  from  a  He-Ne  laser  is  spatially 
filtered  to  create  a  Gaussian  point  source  which  is  then  imaged  in  the  detector  plane 
by  the  lenses.  The  PCGH  is  kinematically  mounted  [3-11]  on  a  micropositioner 
which  affords  three  degrees  of  translation  and  one  of  rotation.  Kinematic  mounting 
allows  removal  for  trimming  and  subsequent  replacement  of  the  PCGH  without  dis¬ 
turbing  the  alignment.  In  the  detector  plane,  the  single  UDT-455  photodetector  is 
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Figure  3-18.  Experimental  Setup  for  Demonstration  and  Evaluation  of  PCGHs.  A 

single  detector  scans  across  the  multiple  output  locations  while  a 
second  detector  monitors  the  undiffracted  beam.  Later 
experiments  substitute  a  single-mode  optical  fiber  for  the  spatial 
filter. 

mounted  on  a  motorized  translation  stage.  The  full  width  of  the  main  lobe  from  a 
minimum-size  facet  in  our  system  :s  2  mm.  One  of  several  square  apertures,  selec¬ 
table  in  size  from  1  to  2  mm,  is  placed  in  front  of  the  detector  to  define  and  limit  its 
effective  area.  A  second,  identical  photodiode  monitors  the  DC  beam  to  compensate 
for  laser  power  fluctuations. 

Alignment  of  the  PCGH  is  straightforward.  With  the  DC  detector  removed,  the  mov¬ 
able  detector  is  positioned  to  receive  maximum  signal  from  the  DC  beam.  The  stage 
is  then  translated,  and  the  PCGH  rotated  to  maximize  one  of  the  other  outputs.  This 
aligns  the  outputs  with  the  path  of  the  detector  and  this  alignment  is  unaffected  by 
subsequent  translation  of  the  PCGH.  Centering  the  PCGH  on  the  optic  axis  is 
achieved  by  equalizing  the  signals  from  four  open  apertures  located  just  outside  the 
hologram  boundaries.  Translation  of  the  hologram  along  the  optic  axis  effectively 
determines  beam  diameter,  since  the  beam  is  slightly  converging  in  this  region. 
Beam  diameter  is  adjusted  to  yield  best  agreement  with  the  design  parameters. 

The  setup  as  described  above  wa?  ’sed  to  evaluate  the  PCGHs  of  Figures  3-14  and 
3-17.  These  holograms  were  each  designed  to  produce  10  equal  intensity  outputs 
when  uniformly  illuminated.  Experimentally,  the  outputs  differed  by  as  much  as  a 
factor  of  two  in  intensity.  Much  of  this  variation  was  attributable  to  non-uniform  il¬ 
lumination  (an  insufficiently  expanded  Gaussian).  We  therefore  proceeded  to  experi-. 
ments  that  would  utilize  the  full  Gaussian  beam.  Out  of  several  PCGHs  designed  for 
Gaussian  illumination,  two  were  ultimately  fabricated. 

Figure  3-19  shows  the  partitioning  by  subfacets  of  a  PCGH  designed  for  Gaussian  il¬ 
lumination  and  10  equal  outputs.  Experimentally,  we  found  the  outputs  to  be  equal 
to  within  20  dB  (Figures  3-20  and  3-21),  a  result  we  consider  quite  good  for  an  optical 
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Figure  3-19.  PCGH  Designed  for  Gaussain  Illumination  (diminishing  to  1/e4  at 
the  comers)  and  10  Equal  Intensity  Outputs.  Symmetry  provides 
some  immunity  to  beam  wander. 


beam,  all  10  first-order  beams,  and  several  second-order  beams 
are  visible.  Note  that  bright  sidelobes  occur  away  from  adjacent 
channels. 
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computing  scheme.  A  discouraging  fact  was  that  day-to-day  stability  of  the  setup 
was  not  much  better  than  the  20  dB.  The  dominant  problem  here  was  beam  wander 
in  the  laser,  which  manifested  itself  as  time  dependent  pointing  of  the  Gaussian 
point  source.  With  this  variability,  improving  the  PCGH  through  trimming  proved 
impossible. 

Our  objective  of  pushing  the  numerical  accuracy  beyond  20  dB  demanded  a  more 
stable  source  of  illumination.  For  this  reason,  we  modified  the  experiment  to  use  a 
single-mode  optical  fiber  as  the  point  source.  A  fiber  offers  distinct  advantages  over 
a  pinhole  spatial  filter;  the  shape  and  pointing  of  the  output  beam  are  unrelated  to 
what  goes  on  at  the  input  end.  Unfortunately,  the  fiber  output  distribution  was  only 
approximately  Gaussian.  There  were  also  difficulties  in  adequately  stripping  clad¬ 
ding  modes  and  in  increased  sensitivity  of  intensity  to  input  alignment. 

When  we  evaluated  the  equal  intensity  output,  Gaussian  illumination  PCGH  of  Fig¬ 
ure  3-19  using  the  fiber  optic  source,  the  accuracy  dropped  from  20  dB  to  about  7  dB 
due  to  the  different  illumination  profile.  The  PCGH  of  Figure  3-22,  designed  for 
Gaussian  illumination  and  nine  non-zero  outputs  spanning  a  40  dB  numerical 
range,  had  a  measured  numerical  range  with  the  fiber  optic  source  of  slightly  under 
37  dB.  This  3  dB  discrepancy  was  due  to  the  two  dimmest  channels  each  occupying  a 
full  subfacet.  Numerical  accuracies  of  the  other  eight  channels  were  about  7  dB.  Al¬ 
though  final  accuracies  were  less  than  what  we  initially  had  hoped  for,  further  re¬ 
finement  of  this  method  could  lead  to  a  vector-matrix  multiplication  scheme  with 
high  accuracy,  large  numerical  range,  and  high  optical  efficiency,  particularly  if  im¬ 
plemented  via  phase  gratings. 
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Figure  3-22.  PCGH  Design  for  Gaussian  Illumination  and  9  Non-Zero  Outputs 

Spanning  a  40  dB  Numerical  Range  (y.  =  0,  ym  =  10m/2).  Channel 
10  occupies  most  of  the  bright  central  portion,  whereas  the 
dimmest  channel  is  located  in  the  corner  where  intensity  is 
diminished  by  1/e4. 
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Section  4 

Fabrication  Techniques  for  increased  Diffraction  Efficiency 


Because  of  their  low  weight,  small  size  and  potentially  easy  replicability,  diffractive 
optical  elements  are  an  attractive  alternative  to  conventional  optical  elements  in 
systems  employing  monochromatic  light.  In  the  two  preceding  sections,  we  have 
shown  that  e-beam  computer  generated  holograms  can  have  very  large  space- 
bandwidth  products  and  produce  arbitrarily  prescribed  wavefronts.  The  holograms 
we  generated  for  optical  testing  and  optical  computing  were  binary,  absorption  type 
(chrome-on-glass)  holograms  with  diffraction  efficiencies  limited  to  a  theoretical 
maximum  of  10.1  percent.  In  this  section,  we  consider  e-beam  fabrication  techniques 
for  producing  computer  generated  holograms  offering  much  higher  diffraction  effi¬ 
ciency. 

To  exhibit  high  diffraction  efficiency  a  hologram  must  utilize  phase  modulation 
rather  than  absorption  modulation.  Phase  modulated  holograms  are  typically  char¬ 
acterized  as  either  thick  or  thin,  although  the  transition  is  a  continuous  one  and  in¬ 
termediate  types  are  common  [4-1].  Thin  phase  holograms  can  oe  treated  as  two- 
dimensional  phase-retardation  plates,  while  thick  phase  holograms  are  three- 
dimensional  recordings  of  interference  patterns. 

Thick  phase  holograms  depend  on  periodic  coupling  along  the  propagation  direction 
such  that  proper  phase  matching  is  achieved  for  only  one  diffracted  order.  When  the 
refractive  index  modulation  and  film  thickness  are  properly  chosen,  nearly  all 
energy  is  coupled  to  this  one  order.  Thick  phase  holograms,  when  interferometrical- 
ly  recorded  in  dichromated  gelatin,  can  exhibit  diffraction  efficiencies  approaching 
100  percent  [4-2].  Unfortunately,  interferometric  recording  lacks  flexibility  for 
producing  arbitrary  wavefronts.  It  is  also  difficult  to  construct  thick  phase  holo¬ 
grams  for  use  at  long  wavelengths  since  efficient  recording  materials  are  not  red  or 
infrared  sensitive. 

One  of  our  objectives  has  been  to  develop  fabrication  techniques  for  producing  high 
diffraction  efficiency  e-beam  computer  generated  holograms.  We  have  restricted  our 
consideration  to  holograms  of  the  surface-relief  type;  while  e-beam  direct  writing  of 
phase  gratings  in  amorphous  chalcogenite  films  has  been  demonstrated,  the  index 
modulation  (An/n~0.01)  and  film  thickness  (<5  micron)  are  insufficient  to  permit 
diffraction  efficiencies  greater  than  a  few  percent.  We  have  investigated  two  types  of 
surface-relief  holograms.  The  first  type  is  a  transmission  hologram  produced  by  ion 
milling  of  a  dielectric  substrate.  The  second  type  is  a  blazed,  reflection  hologram 
produced  by  anisotropic  etching  of  a  silicon  substrate.  The  remainder  of  this  section 
reports  details  of  this  work. 
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ION-MILLED  HOLOGRAPHIC  OPTICS 

With  slight  modification,  our  e-beam  fabrication  procedure  for  chrome-on-glass 
absorption  type  holograms  can  also  produce  thin  phase  holograms  in  e-beam  resist. 
The  modification  involves  placing  the  metallization  (required  for  conductivity  dur¬ 
ing  exposure)  on  top  of  the  e-beam  resist  and  then  removing  it  by  etching  before  re¬ 
sist  development.  The  result  is  a  square-wave  grating  in  e-beam  resist.  Recently,  it 
has  been  demonstrated  that  this  technique  is  also  capable  of  producing  high-efficien- 
cy  (60  to  70  percent)  sawtooth  grating  profiles  if  the  electron  dose  is  suitably  varied 
[4-4].  A  disadvantage  is  that  the  resulting  holograms  are  hygroscopic  and  lack  en¬ 
vironmental  durability. 

Rather  than  use  the  e-beam  resist  directly  as  a  hologram,  we  use  it  as  a  mask  for  ion 
milling  the  grating  pattern  into  the  glass  (or  quartz)  substrate.  This  produces  an  all¬ 
glass  hologram  that  is  very  environmentally  durable.  It  is  also  nonabsorptive,  mak¬ 
ing  it  capable  of  high  diffraction  efficiency  and  a  high  damage  threshold.  The  follow¬ 
ing  paragraphs  discuss  these  ion-milled  holograms  in  detail. 

Theory 

Dielectric  surface-relief  gratings  have  received  extensive  theoretical  treatment  [4-5]. 
At  low  spatial  frequencies  (grating  period  >  >  1  wavelength)  these  holograms  be¬ 
have  as  thin  phase  gratings.  They  typically  produce  many  diffracted  orders  (Figure 
4-1)  whose  intensities  depend  on  the  groove  profile.  A  square-wave  profile,  as 
produced  by  ion  milling,  yields  a  maximum  diffraction  efficiency  of  40.5  percent. 
This  maximum  occurs  for  a  milling  depth  corresponding  to  a  phase  retardation  of  t 
radians  and  a  groove  width  equal  to  half  the  grating  period. 


INCIDENT 

WAVE 


Figure  4-1.  Diffraction  by  a  Low  Spatial  Frequency  Ion-Milled  Grating.  Many 
diffracted  orders  compete  for  incident  light,  with  no  single  order 
receiving  more  than  40.5  percent  of  the  total.  Diffraction  efficiency 
depends  only  weakly  on  angle  of  incidence. 
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DIFFRACTION  EFFICIENCY 


GRATING  DEPTH  (MICRONS) 

Figure  4-2.  Diffraction  Efficiency  Versus  Groove  Depth  for  Ion-Milled 

Holographic  Optics.  A  refractive  index  of  1.5  and  normal  incidence 
have  been  assumed.  Maximum  diffraction  efficiency  occurs  for  a 
milling  depth  of  one  wavelength,  corresponding  to  a  phase  depth  of 
x  radians. 
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DUTY  CYCLE 

Figure  4-3.  Diffraction  Efficiency  Versus  Groove  Width  for  Ion-Milled 

Holographic  Optics.  Maximum  first-order  diffraction  efficiency 
occurs  for  a  groove  width-to-spacing  ratio  of  0.5.  Optimum  groove 
width  is  verified  by  the  absence  of  even  diffracted  orders. 

Figure  4-3  shows  the  dependence  of  diffraction  efficiency  on  groove  width  for  the 
same  grating  parameters  as  Figure  4-2.  The  milling  depth  has  been  optimized  for 
maximum  diffraction  efficiency.  It  is  seen  that  groove  width  must  be  controlled  to  a 
fraction  (- 1/10)  of  the  grating  period  in  order  to  maximize  first-order  diffraction  effi¬ 
ciency.  Absence  of  the  second  diffracted  order  provides  a  sensitive  indicator  of  op¬ 
timum  groove  width.  Figure  4-3  is  invalid  for  either  very  large  or  very  small  groove 
widths  since  the  sine  series  used  to  represent  the  square-wave  groove  profile  was 
truncated  at  ten  terms. 

At  low  spatial  frequencies,  surface-relief  gratings  produce  many  diffracted  orders  as 
was  shown  in  Figure  4-1.  As  grating  period  is  decreased,  fewer  diffracted  orders  are 
available  to  compete  for  the  incident  light.  Also,  as  the  diffraction  angles  become 
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larger,  the  three-dimensional  nature  of  the  groove  profile  becomes  important  and 
these  gratings  can  no  longer  be  treated  as  simple  phase-retardation  plates.  The 
polarization  state  of  the  incident  beam  also  becomes  significant. 

We  expect  maximum  attainable  diffraction  efficiency  to  increase  with  grating 
spatial  frequency.  At  ultra-high  spatial  frequencies,  the  geometry  is  such  that  only 
one  forward-diffracted  order  can  exist  to  compete  for  incident  light;  all  other  orders 
being  evanescent.  This  situation  is  depicted  in  Figure  4-4.  The  minimum  spatial  fre¬ 
quency  to  just  satisfy  this  condition  corresponds  to  a  grating  period  of  1.5 
wavelengths  and  an  incident  angle  d0  =  sin'1  (1/3)  =  19.5  degrees.  Unfortunately, 
our  model  suffers  numerical  inaccuracies  when  groove  depth  becomes  comparable  to 
groove  spacing.  We  have  therefore  been  unable  to  analyze  holograms  of  ultra-high 
spatial  frequency.  However,  there  are  a  number  of  theoretical  and  experimental  re¬ 
sults  to  indicate  that  maximum  diffraction  efficiency  increases  with  increasing 
spatial  frequency.  Moharam  and  Gaylord  [4-5]  have  calculated  a  diffraction  efficien¬ 
cy  of  88.5  percent  for  a  square-wave  grating  of  period  1  wavelength  and  depth  1.55 
wavelengths.  Enger  and  Case  [4-2]  have  demonstrated  88-percent  diffraction  efficien¬ 
cy  for  photoresist  gratings  of  period  0.33  micron  and  depth  0.37  micron.  We  there¬ 
fore  expect  an  increase  in  diffraction  efficiency  as  groove  depth  becomes  comparable 
to  groove  spacing. 


t 


ORATING 


Figure  4-4.  Diffraction  by  an  Ultra-High  Frequency  Ion-Milled  Grating.  The 

geometry  is  such  that  only  one  diffracted  order  is  transmitted;  all 
others  are  evanescent. 
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Design 


Figure  4-5  illustrates  our  basic  design  for  an  ion-milled  holographic  lens.  We  wished 
to  demonstrate  a  high  diffraction  efficiency,  large  numerical  aperture  (NA)  lens 
which  would  focus  a  collimated  He-Ne  laser  beam  to  a  diffraction  limited  point.  The 
diffraction  limited  focal  spot  size  is  Ax  >  X/NA,  where  NA  *  sin  [tan'1  (R/f)].  Fabri¬ 
cation  difficulty  is  related  to  the  smallest  grating  period  Amin.  We  fix  the  maximum 
grating  period  at  Amax  =  2  Amin  in  order  to  assure  angular  separation  of  the  first  dif¬ 
fracted  order,  limit  the  range  of  feature  sizes  to  be  drawn,  and  avoid  possible  lower 
diffraction  efficiency  at  lower  spatial  frequencies.  These  considerations  led  us  to  the 
geometry  depicted  in  Figure  4-5. 


FOCUS 


Figure  4*5.  Design  for  Ion-Milled  Holographic  Lens.  The  lens  focuses  a 
collimated  He-Ne  laser  beam  to  a  point. 


Table  4-1  lists  the  parameters  of  several  holographic  lenses  we  have  designed  along 
the  lines  of  Figure  4-5.  The  quantity  of  data  required  to  specify  one  of  these  lenses  is 
related  to  both  Amin  and  to  hologram  diameter.  We  arbitrarily  selected  a  diameter  of 
5  mm  for  these  lenses.  In  cases  where  the  data  quantities  became  excessive  (all  but 
CGH100  and  CGH101),  we  specified  cylindrical  lenses  with  focal  power  in  only  one 
axis.  Design  CGH110C,  requiring  a  0.5-micron  line  width,  taxed  the  resolution  limit 
of  our  e-beam  lithography  system. 

Fabrication  Techniques 

Figure  4-6  illustrates  the  basic  fabrication  process  for  ion-milled  holographic  lenses. 
The  initial  lithography  is  essentially  the  same  as  for  our  standard  chrome-on-glass 
fabrication  process  depicted  in  Figure  1-3.  A  metallization  layer  is  still  required  for 
electrical  conductivity  during  exposure,  so  we  generally  use  commercially  available 
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Table  4-1.  Ion-Milled  Holographic  Lens  Designs 


Identifier 

Focal  Length 
(mm) 

f-number 

Amin 

(microns) 

Fringes 

CGH100 

160 

f/32 

10.2 

370 

CGH101 

80 

f/16 

5.1 

739 

CGH102 

40 

f/8 

2.5 

1474 

CGH103 

30 

f/6 

1.9 

1959 

CGH108 

25 

f/5 

1.6 

2329 

CGH109 

20 

m 

1.3 

2906 

CGH110 

15 

m 

1.0 

3733 

plates,  which  come  with  a  layer  of  chrome  ( -  650A)  and  resist  (6000APMMA)  al¬ 
ready  deposited.  After  exposure  and  development,  the  resist  acts  as  an  etch  mask  for 
ion  milling  of  the  grating  pattern  into  the  glass  substrate. 

The  relative  milling  rates  of  the  various  materials  determine  the  depth  to  which  we 
can  mill  into  the  substrate.  Milling  is  done  with  argon  ions  in  an  Ion  Tech  dual 
beam  ion  mill.  Relative  milling  rates  of  key  materials  are  indicated  in  Table  4-2. 
Our  standard  glass-chrome-PMMA  mask  plates  can  only  be  milled  to  a  depth  of 
22(M)A  before  the  PMMA  is  entirely  depleted.  From  Figure  4-2,  we  see  that  a  depth 
of  about  6000A  is  required  for  maximum  diffraction  efficiency.  Chemical  etching  of 
the  chrome  from  exposed  regions  would  increase  milling  depth  to  only  2800A 
before  both  PMMA  and  chrome  were  depleted.  Instead,  our  procedure  is  to  first  etch 
through  the  chrome  to  create  a  chrome-on-glass  hologram,  but  then  use  it  as  a  mask 
for  contact  printing  the  hologram  pattern  onto  a  second  glass  plate  coated  with 
7 000-8000 A  of  photoresist.  The  lower  ion  etch  rate  of  photoresist  then  allows  mil¬ 
ling  to  the  required  depth. 

While  the  contact  printing  method  works  well  for  producing  low  spatial  frequency 
holograms,  it  does  not  allow  the  full  resolution  capability  of  e-beam  lithography  to 
be  utilized.  We  have,  therefore,  developed  two  alternative  fabrication  processes  for 
holograms  of  high  spatial  frequency. 

Our  first  alternative  process  replaces  the  PMMA  with  a  Honeywell-proprietary  e- 
beam  resist,  EP25H.  This  resist  has  half  the  ion  etch  rate  and  twice  the  sensitivity 
of  PMMA.  A  6500A  layer  of  EP25H  allows  milling  to  the  optimum  depth  before 
both  resist  and  chrome  are  depleted.  With  faster  etching  quartz  substrates,  the 
thickness  of  EP25H  can  be  reduced  to  about  3000 A,  allowing  still  higher  resolu¬ 
tion  patterns  to  be  drawn. 
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GLASS  SUBSTRATE 


COAT  WITH  RESIST 


EXPOSE  AND  DEVELOP 


ION  MILL 


REMOVE  RESIST 


Figure  4-6.  Fabrication  Process  for  Ion-Milled,  Holographic  Optics.  E-beam 
lithography  is  used  to  produce  a  hologram  in  resist  on  a  glass 
substrate.  The  resist  then  serves  as  a  mask  for  ion  milling  the 
hologram  into  the  glass.  After  removal  of  the  resist,  we  are  left  with 
a  high-efficiency,  all-glass  hologram. 
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Table  4-2.  Ion  Beam  Etch  Rates  (500  eV  Argon  Ions  at  1  mA/cm2,  normal 
incidence) 


Glass  (soda-lime) 

200  A“/min 

Quartz 

400 

PMMA  e-beam  resist 

550 

EP25H  e-beam  resist 

200 

Chrome 

530 

AZ1 350  photoresist 

230 

Our  second  alternative  fabrication  process,  shown  diagrammatically  in  Figure  4-7,  is 
the  ultimate  in  contact  printing.  A  chrome-on-glass  hologram  is  coated  directly  with 
photoresist  and  then  exposed  by  illuminating  from  the  back  side.  The  chrome  serves 
as  an  exposure  mask,  while  the  developed  photoresist  serves  as  an  etch  mask  for  ion 
milling.  As  long  as  the  chrome  remains  intact,  the  process  can  be  repeated  as  many 
times  as  necessary  for  milling  to  any  desired  depth.  This  process  allows  precise 
control  of  milling  depth  since  the  resist  can  be  removed  to  check  groove  depth  and 
then  replaced  for  final  milling  to  desired  depth.  With  quartz  substrates  and  UV 
exposure,  this  technique  can  be  extended  to  very  high  resolution  patterns. 

Experimental  Results 

We  have  been  quite  successful  in  fabricating  low  spatial  frequency  holographic 
lenses  with  diffraction  efficiencies  near  the  theoretical  limit  of  40.5  percent.  Achiev¬ 
ing  a  groove  width-to-spacing  of  0.5  has  proven  to  be  a  greater  challenge  than 
achieving  the  correct  milling  depth,  particularly  for  higher  spatial  frequency  pat¬ 
terns.  This  difficulty  in  controlling  groove  width  prevented  us  from  attafhing  the 
higher  diffraction  efficiencies  we  had  hoped  for  at  grating  periods  near  1  micron. 

Figure  4-8  is  an  SEM  photo  of  the  ion-milled  holographic  lens  CGH100.  This  f/32 
lens  has  a  near  optimum  square-wave  groove  profile,  as  can  be  seen  by  examining 
its  diffraction  pattern  in  Figure  4-9.  The  incident  laser  beam  was  somewhat  larger 
in  diameter  than  the  hologram,  causing  the  bright  ring  around  the  position  of  the 
undiffracted  N  =  0  order.  The  bright  spot  to  the  right  of  this  ring  is  the  focused  N  = 
+1  order  beam,  while  the  bright  disk  to  the  left  is  the  N  =  1  order  beam  for  which 
the  hologram  acts  as  a  negative  lens.  The  fainter  images  to  either  side  are  the  N  = 
±3  order  beams.  The  darkness  of  the  N  =  0  order  and  the  absence  of  any  even  (N  = 
±  2)  orders  attest  to  the  optimum  groove  profile. 
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CHROME-ON-GLASS  HOLOGRAM 


l 


J 


COAT  WITH  PHOTORESIST 


EXPOSE  FROM  BELOW  AND  DEVELOP 


f  tfttmtm 

UV  RADIATION 


ARGON  IONS 

uuimm 

ION  MILL 


Figure  4-7.  Alternative  Fabrication  Process  for  Ion-Milled  Holographic  Optics. 

This  process  allows  deeper  milling  since  the  photoresist,  which  has 
a  lower  milling  rate  than  e-beam  resist,  may  be  applied  repeatedly. 


#0‘ 


Figure  4-9.  Diffraction  Pattern  from  Ion-Milled  Holographic  Lens  CGH100.  The 

center  of  the  ring  is  the  N  =  0  order,  while  the  bright  spot  to  the 
right  is  the  focused  N  =  +  1  order. 


Figure  4-10  is  an  SEM  photo  of  the  cylindrical  holographic  lens  CGH110C.  This 
lens,  with  a  grating  period  of  1.0  micron,  is  the  highest  spatial  frequency  lens  we 
have  fabricated.  Control  of  linewidth  in  the  relatively  thick  (5000A)  resist  during 
exposure  and  subsequent  ion  milling  proved  exceedingly  difficult.  The  small  feature 
sizes  required  that  each  sample  be  produced  by  e-beam  direct-writing.  This  limited 
the  number  available  for  process  development.  Of  several  samples  fabricated,  none 
had  a  diffraction  efficiency  exceeding  20  percent. 

To  summarize,  we  found  that  we  could  fabricate  ion-milled  holographic  lenses  hav¬ 
ing  diffraction  efficiencies  near  the  theoretical  maximum  of  40.5  percent  if  grating 
periods  were  greater  than  about  2.5  microns.  Theory  predicts  an  increase  in  diffrac¬ 
tion  efficiency  for  grating  periods  of  1  micron  and  under,  but  difficulties  in  control¬ 
ling  the  grating  profile  at  such  small  geometries  resulted  instead  in  reduced 
diffraction  efficiency. 

BLAZED,  CRYSTALLINE  HOLOGRAPHIC  OPTICS 


Our  second  approach  to  achieving  high  diffraction  efficiency  holographic  optical  ele¬ 
ments  uses  anisotropic  etching  of  crystals  to  produce  blazed,  reflection  holograms. 
Blazed  diffraction  gratings  are  often  made  by  ruling  techniques  and  can  have  dif¬ 
fraction  efficiencies  approaching  100  percent.  Anisotropic  etching  techniques  are 
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Figure  4-10.  SEM  Photo  of  Ion-Milled  Holographic  Lens  CGH110C.  The  groove 
spacing  is  about  1  micron. 

used  to  create  grating  patterns  for  distributed  feedback  laser  diodes  [4-8].  Figure 
4-11  shows  one  such  grating.  Our  approach  has  been  to  combine  e-beam  lithography 
with  anisotropic  etching  to  produce  blazed,  holographic  lenses  in  silicon  substrates. 

Theory  and  Fabrication  Techniques 

The  anisotropic  etching  of  silicon  is  a  technique  of  growing  importance  for  the  fabri¬ 
cation  of  both  integrated  circuits  [4-9]  and  micromechanical  structures  [4-10,  4-11]. 
Several  ternary  liquid  etchants  have  been  reported  that  etch  silicon  at  different 
rates  in  the  principle  crystallographic  directions.  These  consist  of  an  oxidant  such  as 
potassium  hydroxide  (KOH),  which  oxidizes  silicon  to  hydrated  silica,  and  a  complex- 
ing  agent  such  as  iso-2-propyl  alcohol  (IPA),  which  reacts  with  the  silica  to  form  a 
soluble  complex  ion  [4-12].  Figure  4-12  indicates  the  principle  crystallographic 
planes  of  (100)  silicon.  The  (111)  planes  are  the  slowest  etching  planes.  The  (110)  and 
(100)  planes  can  etch  400  times  faster  than  the  (111)  planes.  Whether  the  (110)  plane 
etches  faster  than  the  (100)  plane,  or  vice  versa,  depends  upon  the  etch  formulation 
used.  The  fastest  etching  plane  is  the  (211)  plane.  The  etchant  does  not  attack 
or  Si02. 

The  fabrication  process  for  blazed,  crystalline  grating  structures  is  depicted  in  Fig¬ 
ure  4-13.  E-beam  lithography  is  used  to  delineate  a  grating  pattern  in  a  Si3N4  or 
Si02  etch  mask  on  (100)  silicon.  The  grooves  must  be  nominally  parallel  to  the  slow 
etching  (111)  planes.  Etching  proceeds  in  the  (100)  direction  (down)  and  is  stopped  by 


Figure  4-11.  SEM  Photo  of  Anisotropically  Etched  Grating  in  Gallium  Arsenide. 

The  cleaved  edge  view  shows  the  vee-groove  profile.  This  0.376- 
micron  linear  grating  was  interferometrically  recorded. 


64.7  DEG 


Figure  4-12.  Orientation  of  Major  Crystallographic  Planes  in  (100)  Silicon 
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Si  OR  GaAs  SUBSTRATE 


COAT  WITH  Si3N4  AND  RESIST 


EXPOSE  AND  DEVELOP 
ETCH  Si3N4 


PREFERENTIAL  Si  OR 
GaAs  ETCH 

ANISOTROPIC  ETCHING  OF 
Si  OR  GaAs 


REMOVE  Si3N4 
AND  METALLIZE 


Figure  4-13.  Fabrication  Process  for  Blazed,  Crystalline  Holographic  Optics.  E- 
beam  lithography  is  used  to  produce  a  hologram  in  Si3N*  on  (100) 
silicon.  The  narrow  lines  of  resist  must  run  roughly  parallel  to  the 
(111)  planes.  The  anisotropic  etch  is  stopped  by  the  Si3NU  mask 
and  the  (111)  planes,  leaving  precise  vee-grooves.  After  cleaning, 
the  hologram  is  metallized  or  used  as  a  master  for  embossing 
copies. 
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the  mask  and  the  (111)  planes,  leaving  precise  vee-grooves  with  faces  inclined  at 
54.7  degrees  to  the  substrate  normal. 

Diffraction  by  a  blazed,  crystalline  holographic  grating  is  depicted  in  Figure  4-14. 
Diffraction  angles  are  determined  by  the  groove  spacing  and  the  grating  equation. 
Tilted  groove  faces  preferentially  reflect  light  into  particular  diffracted  orders. 
When  both  the  grating  equation  and  the  reflection  law  can  be  satisfied  simultane¬ 
ously,  much  of  the  incident  energy  is  channeled  into  a  single  diffracted  order. 

As  blaze  angles  go,  the  54.7-degree  blaze  produced  by  etching  of  (100)  crystals  is  very 
steep.  Other  blaze  angles  are  possible,  but  would  require  specially  cut  crystals.  Very 
high  blaze-angle  gratings,  because  of  their  symmetric  groove  profiles,  show  peculiar 
first-order  diffraction  efficiencies  [4-13].  Diffraction  efficiency  is  found  to  depend  sig¬ 
nificantly  on  polarization  state  and  on  departure  from  Littrow  configuration  (zero 
angular  deviation). 

The  sharply  peaked  grating  profile,  as  depicted  in  Figure  4-14,  poses  a  number  of 
fabrication  difficulties.  The  linewidth  of  the  etch  mask  pattern  must,  of  course,  be 
sufficient  that  undercutting  does  not  free  the  mask  before  the  vee-grooves  are  etched 
to  their  full  depth.  Any  excess  width  will  result  in  a  flat-topped  grating  profile  which 
reduces  diffraction  efficiency  by  preferentially  directing  light  into  the  zero-order  (un¬ 
diffracted)  beam.  Flat  tops  become  flat  bottoms  if  the  etched  hologram  is  used  as  a 
master  for  embossing  copies.  Flat  bottoms  filling  less  than  1/3  of  the  groove  depth 
will  be  shadowed  by  adjacent  peaks. 

Design 

There  are  a  number  of  design  constraints  associated  with  blazed,  crystalline  holo¬ 
graphic  optics.  These  can  be  attributed  primarily  to  the  fixed  blaze  angle  and  groove 
orientation  imposed  by  the  crystal  geometry. 

Figure  4-15  shows  our  designs  for  blazed,  crystalline  holographic  lenses.  The  central 
ray  is  oriented  so  that  it  will  be  specularly  reflected  off  the  groove  faces.  In  addition, 
we  require  that  the  incident  angle  be  shallower  than  the  diffraction  angle.  This  pre¬ 
vents  rays  from  striking  the  opposite  face  of  a  groove  as  they  leave  the  grating  and 
provides  our  rationale  for  distinguishing  between  collimators  and  condensers.  To 
maximize  numerical  aperture,  we  make  the  collimated  beam  as  large  as  possible 
without  engulfing  the  point  source  or  focus. 

Blazed  gratings  can  be  designed  to  operate  in  nearly  any  diffractive  order.  Very  high 
blaze-angle  gratings  are  rarely  used  in  first  order  [4-13].  The  advantage  of  a  higher 
order  is  that  the  grating  period  can  be  larger  (although  the  required  spatial  accuracy 
is  unaffected).  A  disadvantage  is  that  adjacent  orders  are  closer  in  angle  and  may 
compete  for  the  light.  This  is  a  problem  particularly  for  our  blazed,  crystalline  holo¬ 
graphic  optics,  since  the  grating  period  will  vary  but  the  blaze  angle  will  not.  There¬ 
fore,  we  might  expect  the  largest  numerical  aperture  when  working  in  the  first  or¬ 
der.  However,  the  small  grating  periods  required  (X/V  6  in  Littrow  configuration) 
make  first-order  blazed,  crystalline  gratings  impractical. 
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into  a  particular  diffracted  order. 
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Figure  4-16.  Two  Designs  for  Blazed,  Crystalline  Holographic  Optics.  The 
geometry  has  been  chosen  so  that  the  central  ray  will  be 
specularly  reflected  off  the  groove  faces,  which  are  inclined  at 
54.7  degrees.  In  order  that  light  not  undergo  multiple  reflections, 
the  incident  angle  must  be  steeper  than  the  diffraction  angle.  This 
criterion  fixed  the  point  source  (or  focus)  at  the  edge  of  the 
collimated  beam. 
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Table  4-3  lists  parameters  of  several  blazed,  crystalline  holographic  lens  designs. 
These  holograms  are  5  mm  in  diameter  and  designed  for  operation  at  the  He-Ne 
(6328 A.)  wavelength.  In  each  case,  the  diffractive  order  was  chosen  as  large  as 
possible  without  having  the  blaze  angle  favor  adjacent  orders  over  some  portion  of 
the  hologram.  A  line  width-to-spacing  ratio  of  2/3  was  chosen  for  these  holographic 
lenses  with  the  intention  of  replicating  by  embossing. 


Table  4-3.  Blazed,  Crystalline  Holographic  Lens  Designs 


Identifier 

Focal 

Length 

(nun) 

f-number 

H 

Order 

Type 

Fringes 

CGH104 

80 

f/16 

2.7 

7th 

Condenser 

1841 

CGH105 

80 

f/16 

2.7 

7th 

Collimator 

1836 

CGH106 

40 

f/8 

1.5 

4th 

Condenser 

3382 

CGH107 

40 

08 

1.5 

4th 

Collimator 

3222 

Experimental  Results 

Our  attempts  to  fabricate  blazed,  crystalline  holographic  lenses  were  plagued  by  un¬ 
dercutting  of  the  Si3N4.  While  some  undercutting  was  expected,  the  magnitude  of 
the  effect  was  not  well  understood. 

As  an  example,  we  describe  the  fabrication  of  the  lens  CGH104.  This  lens  was  de¬ 
signed  as  a  7th  order  f/16  condenser  with  a  minimum  grating  period  of  2.7  microns. 
The  pattern  was  drawn  as  a  chrome-on-glass  mask  with  a  line  width-to-spacing  ratio 
of  2/3.  This  hologram  was  then  printed  in  AZ1450B  photoresist  on  (100)  silicon 
substrates  coated  with  a  1500A  layer  of  Plasma  etching  was  used  to  open 

grooves  in  the  S^N*. 

The  anisotropic  etchant  used  for  CGH104  consisted  of  160  gm  KOH  in  110  ml  IPA 
and  200  ml  water.  Figure  4-16  shows  the  sample  after  etching  for  12  minutes  at 
35  °C.  In  Figure  4-  16a,  undercutting  of  the  Si3N4  mask  is  nearly  complete,  and  yet 
the  vee-grooves  still  do  not  have  sharp  bottoms.  Also,  the  sides  of  the  grooves  show 
numerous  ledges.  Figure  4-  16b  shows  another  portion  of  the  hologram  where  the 
mask  has  lifted  off  completely  and  some  of  the  peaks  have  begun  to  erode.  Still,  the 
grooves  do  not  have  sharp  bottoms. 

Some  undercutting  of  the  mask  and  ledge  formation  on  the  (111)  faces  were  expected 
due  to  misorientation  of  the  grooves  with  respect  to  the  true  (111)  planes  [4-14].  For 
a  perfectly  oriented  sample,  with  grooves  parallel  to  the  (111)  planes,  the  rate  of 
downward  etching  is  r(10Q),  while  the  rate  of  undercutting  is  ru  =  >/  (3/2)  r(111).  Here 
r(joo)  and  r(iii)  are  etch  rates  in  the  (100)  and  (111)  directions  respectively.  If  grooves 
are  misoriented,  then  ledges  will  form  and  etching  of  these  ledges  can  contribute  sig¬ 
nificantly  to  undercutting. 
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b)  In  another  portion  of  the  hologram,  the  Si3N4  mask  has  completely  lifted 
off  and  some  of  the  peaks  have  begun  to  erode.  Misorientation  ledges  are 
clearly  visible  on  the  groove  sides. 


Figure  4*16.  Blazed,  Crystalline  Holographic  Lens  CGH104  After  Etching  for  12 
Minutes  in  KOH  Solution 
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For  grooves  etched  in  (110)  silicon,  Kendall  [4-14]  finds  the  etch  rate  of  nominal  (111) 
faces  to  be  r(111)  +  1.6  r(110)  A9  where  A0  is  the  angular  misalignment  and  r(110)  > 
400  rniiv  Applying  Kendall’s  results  to  our  geometry  yields  a  rate  of  undercutting 
ru  =  v  (3/2)  r(111)  +  1.6  r(UO)A0.  The  vee*grooves  should  reach  their  full  depth  before 
undercutting  is  complete  if  r(10o)  >  3V  2  ru.  Adopting  r(l(y),  =  r(U0)  =  400  r(111)  we 
find  that  complete  undercutting  should  not  occur  unless  Ab  <  35  degrees.  Our  holo¬ 
graphic  lenses  have  some  misorientation  of  grooves  with  respect  to  the  true  (111) 
planes  because  of  groove  curvature,  but  this  misorientation  is  only  ±  1.1  degree  for 
the  f/16  lenses  and  ±  2.2  degrees  for  the  f/8  lenses.  We  therefore  conclude  that  the 
undercutting  observed  cannot  be  entirely  attributed  to  misorientation. 

SUMMARY 

We  have  demonstrated  all-glass  computer  generated  holograms  fabricated  by  e-beam 
lilthography  and  ion  milling.  These  holograms  have  square-wave  groove  profiles,  re¬ 
sulting  in  moderately  high  diffraction  efficiencies  of  40  percent.  They  are  very  en¬ 
vironmentally  durable  and  have  a  high  damage  threshold.  They  can  be  designed  to 
produce  wavefronts  of  arbitrary  complexity,  limited  only  by  e-beam  resolution  and 
hologram  size. 

Attempts  to  increase  diffraction  efficiency  to  80  percent  by  reducing  the  grating  pe¬ 
riod  to  near  1  micron  were  unsuccessful  because  of  difficulties  in  controlling  the  nar¬ 
row  linewidths  in  thick  resist  layers.  In  any  event,  the  applicability  of  such  ultra- 
high  spatial  frequency  holograms  is  marginal  because  of  design  constraints  imposed 
by  their  high  angular  selectivity. 

We  also  attempted  to  fabricate  high-efficiency,  blazed  holograms  using  anisotropic 
etching  of  silicon.  Undercutting  of  the  e-beam  delineated  etch  mask  prevented  us 
from  achieving  the  desired  groove  profile.  Reasons  for  this  occurance  are  not  fully 
understood.  Blazed,  crystalline  holograms  also  suffer  severe  design  constraints  be¬ 
cause  of  their  fixed  blaze  angle  and  groove  orientation. 

General  purpose,  high-efficiency  computer  generated  holograms  will  probably  re¬ 
quire  a  sawtooth  grating  profile.  Others  have  had  some  success  in  producing  such 
profiles  in  PMMA,  both  by  e-beam  lithography  [4-4]  and  by  ion  milling  [4-15].  We 
feel  there  is  still  much  potential  for  fabricating  similar  structures  in  such  durable 
materials  as  glass  or  quartz.  At  present,  however,  the  need  is  for  space-bandwidth 
product  rather  than  high  diffraction  efficiency. 
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Section  5 

Summary  of  Accomplishments 


HOLOGRAM  WRITING  CAPABILITY 

•  Achieved  e-beam  writing  with  a  line-centering  accuracy  of  ±0.2  microns,  or 
0.004  percent  over  a  1-cm  diameter  image. 

•  Identified  optical  plotter  distortion  (DICOMED  Image  Recorder)  to  be  0.05  per¬ 
cent,  with  usable  resolution  of  2000  pixels  per  diameter. 

•  Developed  versatile  and  efficient  encoding  algorithms  permitting  routine  fabri¬ 
cation  of  e-beam  computer  generated  holograms. 

APPLICATIONS  OF  E-BEAM  COMPUTER  GENERATED  HOLOGRAMS 

•  Demonstrated  use  of  e-beam  computer  generated  holograms  in  testing  of  4-inch 
f/2  parabola  and  asymmetric  deep  aspheric  surfaces. 

•  Supplied  about  50  e-beam  computer  generated  holograms  to  Honeywell  Electro- 
Optics  Division  for  aspheric  testing  of  diamond-turned  surfaces. 

•  Proposed  and  developed  an  optical  vector-matrix  multiplication  scheme  utili¬ 
zing  spatially  partitioned,  computer  generated  holograms  (PCGHs)  to  achieve 
high  numerical  accuracy  and  large  numerical  range. 

•  Demonstrated  PCGHs  with  20-dB  numerical  accuracy  and  37-dB  numerical 
range  in  1  x  10  matrix  operations. 

HIGH  DIFFRACTION  EFFICIENCY.  COMPUTER  GENERATED  HOLOGRAMS 

•  Developed  fabrication  techniques  for  producing  40-percent  efficient,  en¬ 
vironmentally  durable,  all-glass  holograms  with  grating  periods  as  small  as  2.5 
microns. 

•  Investigated  anisotropic  etching  techniques  for  producing  blazed,  reflection  hol¬ 
ograms  in  silicon. 
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Appendix 

Miscellaneous  Information 


PERSONS  WHO  PARTICIPATED  IN  THE  RESEARCH 

S.M.  Arnold,  Principal  Investigator 

K.M.  Leung,  Principal  Investigator  (Jan  1980  to  May  1981) 

F.M.  Schmit,  Senior  Research  Scientist 
J.C.  Lindquist,  Research  Scientist 

H.M.  Zimmermann,  Research  Scientist,  Honeywell  Solid  State  Electronics  Division 

S.K.  Case,  Consultant,  Assistant  Professor  of  Electrical  Engineering  at  the  Univer¬ 
sity  of  Minnesota 

J. C.  Wyant,  Consultant,  Professor  of  Optical  Sciences  at  the  University  of  Arizona 

PUBLICATIONS 

K. M.  Leung,  J.C.  Lindquist  and  L.T.  Shepherd,  “E-beam  Computer  Generated  Holo¬ 
grams  for  Aspheric  Testing,”  Proc.  Soc.  Photo-Opt.  Instrum.  Eng.  215,  70-75,  (1980). 

S.M.  Arnold  and  S.K.  Case,  “E-beam  Generated  Holographic  Masks  for  Optical  Vec¬ 
tor-Matrix  Multiplication,”  presented  at  NASA  conference  on  Optical  Information 
Processing  for  Aerospace  Applications,  NASA  Conf.  Pub.  2207,  309-317  (1981). 

S.M.  Arnold  and  S.K.  Case,  “E-beam  Generated  Holographic  Masks  for  Optical  Vec¬ 
tor-Matrix  Multiplication,”  accepted  for  publication  in  Opt.  Eng. 

OTHER  ACTIVITIES 

Transferred  technology  of  aspheric  testing  using  computer  generated  holograms  to 
Honeywell  Electro-Optics  Division  in  March  of  1980. 
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